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We tabulate the ground-state odd-proton and odd- neutron spins and parities, 
proton and neutron pairing gaps, binding energy, one- and two- neutron sepa- 
ration energies, quantities related to /3-delayed one- and two-neutron emission 
probabilities, /3-decay energy release and half-life with respect to Gamow- Teller 
decay, one- and two-proton separation energies, and a-decay energy release and 
half- life for 8979 nuclei ranging from ^^O to ^^^136 and extending from the proton 
drip line to the neutron drip line. Single-particle level diagrams and other quan- 
tities are also presented in graphical form. The starting point of our present work 
is a study of nuclear ground-state masses and deformations based on the finite- 
range droplet model and folded- Yukawa single-particle potential published in a 
previous issue of Atomic Data and Nuclear Data Tables. The /^-delayed 
neutron-emission probabilities and Gamow- Teller /3-decay rates are obtained from 
a quasi-particle random-phase approximation with single-particle levels and wave 
functions at the calculated nuclear ground-state shapes as input quantities. 



*This paper is dedicated to the memory of our friend and colleague Vilen M. 
Strutinsky, who through his years of devoted research on nuclear-structure 
models, most notably the Strutinsky shell-correction method, made possible 
many of the calculations discussed here. 
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1 INTRODUCTION 

In a previous issue of Atomic Data and Nuclear Data Tables we presented a calculation 
of nuclear ground-state masses and deformations for 8979 nuclei ranging from ^^O to ^'^^136 
and extending from the proton drip line to the neutron drip line. ^) The 1992 version of the 
finite-range droplet model and folded- Yukawa single-particle potential that was the basis for 
this calculation is referred to as the FRDM (1992). We here use these ground-state masses and 
deformations as starting points for calculations of additional ground-state properties that are 
useful for astrophysical and other applications. 

An important feature of a mass model is its reliability for nuclei beyond the region used for 
the determination of the model constants. In particular, can one expect the model to be reliable 
for nuclei very far from /5-stability and in the region of superheavy elements? In our mass 
paper ^) we addressed the model reliability for new regions of nuclei by comparing predictions 
of masses that were not included in the data set to which the model constants were determined 
to new experimental data. Since at that time we had available very few new masses over and 
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above the 1989 data set^) from which the model constants were determined, we had to test 
model reliability by simulation. In one simulation, where the model was adjusted to a 1977 
experimental data set, it was found that 351 new masses measured between 1977 and 1989 
were calculated with an increase of only 2% in the model error. In another simulation, where we 
adjusted the model constants only to nuclei in the region Z, N >28 and A < 208, we found that 
there was no increase in the mean error but some increase in the standard deviation for nuclei 
beyond A = 208 that were not included in the limited adjustment. However, for superheavy 
nuclei such as ^^^110 and ^^"^110, the difference in mass predictions between the FRDM (1992), 
whose constants were determined by including nuclei up to A = 263, and the limited adjustment 
was only of the order of 1 MeV. It is significant that these superheavy nuclei are 80 mass units 
heavier than the heaviest nuclide included in the limited adjustment. 

Because our previous studies indicated a very good model reliability for new regions of 
nuclei, we here present calculations of additional nuclear ground-state properties based on the 
same model and the same values of model constants, for the same set of 8979 nuclei considered 
in our mass calculation. ^) Specifically, we consider the following quantities: 

Odd-nucleon spins and parities: 

Projection of the odd-proton angular momentum along the symmetry axis 

and parity of the wave function Op 

Projection of the odd- neutron angular momentum along the symmetry axis 

and parity of the wave function f^JJ 



Lipkin-Nogami pairing gaps: 

Proton pairing gap ^LNp 

Neutron pairing gap ^LNn 

FRDM mass-related quantity: 

Total binding energy -E^bind 

Neutron separation energies: 

One-neutron separation energy S'ln 

Two-neutron separation energy 5211 



Beta-decay properties: 

Probability for producing a final nucleus with mass number A 

following P decay and delayed neutron emission Pa 
Probability for producing a final nucleus with mass number A — 1 

following /? decay and delayed neutron emission Pa-i 
Probability for producing a final nucleus with mass number A — 2 

following j3 decay and delayed neutron emission Pa-2 

Energy released in /3 decay Qp 

Half-life with respect to Gamow- Teller /3 decay Tp 
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Proton separation energies: 

One-proton separation energy Sip 

Two-proton separation energy S'2p 

Alpha-decay properties: 

Energy released in a decay Qa 

Half-life with respect to a decay Ta 



The details of the calculations are given in Sec. 2. Separation energies and energy releases are 
readily obtained from mass differences. The /3-decay half-lives and /3-delayed neutron-emission 
probabilities are obtained from a quasi-particle random-phase approximation (QRPA). In the 
QRPA the single-particle energies and wave functions at the calculated ground-state deformation 
serve as the starting point. The tabulated results are described in Sec. ^. 

After the submission of our mass paper ^) a new mass evaluation ^'^) has become available. 
It contains 217 new masses that were not included in the 1989 data set ^) from which the values 
of the FRDM (1992) constants were determined. Therefore, we are now able to assess the 
reliability, without simulations, of the FRDM (1992) and of several other models that are also 
commonly used in astrophysical calculations. These reliability issues are discussed in Sec. ^. 

As one particular application, we discuss in Sec. |5| the use of the calculated quantities in 
astrophysical rp- and r-process calculations. However, our results also have wider applicability 
to several other areas of astrophysics and to such fields as reactor physics. Section ^ contains 
calculated proton and neutron single-particle level diagrams for representative spherical and 
deformed nuclei throughout the periodic system. 

2 CALCULATIONAL DETAILS 

The quantities studied in this paper are obtained in four different ways. 

1. The odd-proton spin and parity Jlp, odd-neutron spin and parity JIJ^, proton pairing gap 
AlNp, and neutron pairing gap ALNn are microscopic quantities obtained simultaneously 
with the calculated ground-state masses and deformations. They were not published in 
our mass paper because of space limitations. 

2. The total binding energy ii^bind) one- neutron separation energy 5in, two-neutron separa- 
tion energy 52n, /3-decay energy release Qp^ one-proton separation energy Sip, two-proton 
separation energy S'2p, and a-decay energy release Qa are obtained from appropriate dif- 
ferences of the calculated mass excesses. For convenient access we publish them here. 

3. The /3-delayed occupation probabilities Pa, Pa-i, and Pa-2 and /3-decay half-lives Tp are 
obtained from a microscopic quasi-particle random-phase approximation (QRPA). 

4. The a-decay half-life is obtained from the semi-empirical relationship of Viola and 
Seaborg, ^) with constants determined by Sobiczewski, Patyk, and Cwiok. ^) 

2.1 Odd-nucleon spin and parity 

The odd-nucleon spin is simply the projection of the angular momentum along the symmetry axis 
(0 quantum number) for the last occupied proton or neutron level, when this level is occupied 
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by a single nucleon. For odd-proton or odd-neutron nuclei the spin of the nucleus is simply VL^ 
or r^n, respectively. The superscript vr gives the parity of the wave function. 

For spherical nuclei with degenerate levels, the nuclear spin is defined as the maximum value 
of jzi which is Thus, for a spherical nucleus we cannot use as an odd-even spin assignment 
the O value automatically provided for the last occupied single-particle level, since this level is 
randomly assigned any J7 value in the range 1/2 to For slightly deformed nuclei one could 
in principle use deformed assignments, but in practice it would be unrealistic to list a deformed 
assignment for a nucleus with a calculated deformation of, say, €2 = 0.01. We therefore proceed 
in the following manner. For nuclei with a deformation |e2| > ecrit a deformed assignment is 
used. For nuclei with |e2| < Ccrit we calculate the levels for a spherical shape and adopt the 



spherical spin assignment thus obtained. When we compare below in Sec. 3.1 calculated odd- 
particle spins and parities with experimental data, we show that the results are quite insensitive 
to the exact value of ecrit- We choose here 

ecrit = 0.15 



For nuclei whose ground states are calculated to have the octupole shape parameter ea 7^ 
parity is not conserved and is therefore not tabulated. 

2.2 Pairing gaps 

In an extensive study of nuclear pairing'') we investigated both a macroscopic pairing model 
and a microscopic pairing model, which was solved in both the BCS^^^^) and Lipkin-Nogami 
(LN) 12-14-^ approximations. For each model we determined a preferred form of the effective 
pairing interaction and optimum values of the constants of the effective pairing interaction, 
which were obtained from a least-squares minimization of the difference between calculated 
pairing gaps and experimental odd-even mass differences. 

An important result of our previous study is that it is crucial to differentiate between several 
pairing-gap concepts. The most simple concept is the average pairing gap A, which is an 
algebraic relationship such as c/^/A, where c is a constant and A is the number of nucleons in 
the system being studied. The average pairing gap may be regarded as a macroscopic model for 
the nuclear pairing gap, and it may therefore be directly compared with experimental odd-even 
mass differences. 

When a microscopic approach is used the situation is considerably more complicated. In this 
case the quantities that are compared to experimental odd-even mass differences are obtained as 
solutions to microscopic pairing equations, for example the BCS or LN equations. In the BCS 
method it is A that should be directly compared to the odd-even mass differences. However, in 
the LN approximation it is the sum of the pairing gap A and the number-fluctuation constant 
A2, where A and A2 are obtained as solutions of the LN equations, that should be compared to 
odd-even mass differences. We denote this sum by Aln- Thus, Aln = A + A2. 

To solve the usual pairing equations^) one needs in addition to single-particle energies also 
the value of the pairing-strength constant G. This constant depends in a complicated way on 
the number of levels included in the calculation and on the particular nuclear region consid- 
ered. However, it may be determined from an effective-interaction pairing gap Ac by use of a 
Strutinsky-like procedure. ^) At first sight this may seem an unnecessary complication, but the 
advantage is that A^ does not depend on the particular truncation of the single-particle level 
spectrum that is chosen in the calculation. Furthermore, it depends in a very simple way on Z 
and A'^. Therefore, a significant simplification is achieved if one considers Ac to be the primary 
input quantity for pairing calculations, with the constants that enter the function that defines 
Ag to be the pairing-model effective-interaction constants. 
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In our earlier study ^) we obtained the following preferred functional form for the effective- 
interaction pairing gap Ac-: 



Acn = (1) 

Here Z and N are the numbers of protons and neutrons, respectively, / = {N — Z)/{N + Z) is 
the relative neutron excess, and is the surface area of the nucleus divided by the surface area 
of the spherical shape. From root-mean-square minimizations we obtained results consistent 
with t = for both the BCS and LN models. For these cases Eq. (|l|) simplifies to 



iVl/3 



and the effective-interaction pairing gap Aq is determined by one constant for the entire nuclear 
chart, for both protons and neutrons. 

In our nuclear mass calculation^) we performed a refined determination of the effective 
pairing-interaction constant r, with the result that we adopted the value r = 3.2 MeV instead of 
the earlier value r = 3.3 MeV. For details we refer to these earlier studies. ^''^) Although we 
revised the effective-interaction pairing constant r by 3%, the earlier extensive pairing study ^) 
can still serve as an excellent guide to the properties of our current pairing model. Below we 
present further results obtained in the current model. 



2.3 Total binding energy 

The total binding energy £'bind(-^) ^) is related to the atomic mass excess M(Z, N) through the 
simple relationship 

£;bmd = ZMu + iVMn - M{Z, N) (3) 

where Mh is the hydrogen-atom mass excess and is the neutron mass excess. The total 
binding energy includes the binding energy of the Z electrons comprising the atom, which we 
approximate by a^Z'^-^^, with a^i = 1.433 x 10"^ MeV. 

For the benefit of workers in other fields, who are often confused by the conventions adopted 
in atomic masses, we mention that the reason that the atomic mass excess is tabulated instead 
of the atomic mass itself is simply to eliminate the repetitive tabulation of additional leading 
numbers that can easily be restored by adding the mass number A times the mass unit u, which 
is 1/12 the mass of the ^^C atom, to the tabulated quantity. Also, the reason that the atomic 
mass is considered rather than the nuclear mass is that the former is the actual experimentally 
measured quantity, whereas the latter is less accurate because its extraction requires a knowledge 
of the binding energy of the Z atomic electrons. 

For those applications where it is necessary to know the actual mass of the nucleus itself, 
its value (in MeV) can be found from the atomic mass excess tabulated in Ref. ^) by use of the 
relationship 

Mnucleus = Au + M - ZrUe + UelZ^'^^ (4) 

where mg = 0.51099906 MeV is the mass of the electron, i^'^^) As discussed in Ref. ^), the value 
u = 931.5014 MeV that was used in the interim 1989 mass evaluation^) should be used for the 
atomic mass unit in Eq. 
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2.4 Neutron separation energies 

The one- and two-neutron separation energies Sin{Z,N) and S2niZ,N) are obtained from the 
mass excesses through the differences 

Sin{Z,N) = M{Z,N -1) + M^-M{Z,N) 

S2n{Z,N) = M{Z,N -2) + 2Mn- M{Z,N) (5) 



2.5 P-decay properties 

The formahsm we use to calculate Gamow- Teller (GT) /3-strength functions is fairly lengthy, 
since it involves adding pairing and Gamow- Teller residual interactions to the folded- Yukawa 
single-particle Hamiltonian and solving the resulting Schrodinger equation in the quasi-particle 
random-phase approximation. Because this model has been completely described in two previous 
papers, ^^'^^) we refer to those two publications for a full model specification and for a definition 
of notation used. We restrict the discussion here to an overview of features that are particularly 
relevant to the results discussed in this paper. 

It is well known that wave functions and transition matrix elements are more affected by 
small perturbations to the Hamiltonian than are the eigenvalues. When transition rates are 
calculated it is therefore necessary to add residual interactions to the folded- Yukawa single- 
particle Hamiltonian in addition to the pairing interaction that is included in the mass model. 
Fortunately, the residual interaction may be restricted to a term specific to the particular type 
of decay considered. To obtain reasonably accurate half-lives it is also very important to include 
ground-state deformations. Originally the QRPA formalism was developed for and applied 
only to spherical nuclei, i^'^'') The extension to deformed nuclei, which is necessary in global 
calculations of /3-decay properties, was first described in 1984. ^'^) 

To treat Gamow- Teller f3 decay we therefore add the Gamow- Teller force 

Vgt = 2xGT ■■ P^- • f3^+ ■■ (6) 

to the folded- Yukawa single-particle Hamiltonian, after pairing has already been incorporated, 
with the standard choice XGT = 23 MeV/^. ^''~^*^) Here (3^-^= J2i^iti~ the Gamow- Teller 
/3^-transition operators, and the colans mean that all contractions in the quasi-particle represen- 
tations of the enclosed operator are to be ignored. The correlations generated by the GT force 
are of specific importance to the Gamow- Teller decays, which are the dominant decay modes in 
many nuclei of astrophysical interest. Other types of residual interactions are of importance for 
other decay modes, but leave the Gamow- Teller decay rates unaffected, and can consequently 
be ignored for our present purpose. 

It should be noted that the RPA treatment formulated by Halbleib and Sorensen^^) incor- 
porates only particle-hole correlations of specific importance to GT transitions. It has been 
proposed ^^'^^) that the effect of neglected particle-particle terms may be significant for tran- 
sitions. We later address this question in Sec. Moreover, the RPA treatment may not contain 
enough ground-state correlations.^^) However, in view of the present uncertainties regarding 
these points we leave possible further refinements for future consideration. Some additional 
comments are made in Sec. |3.5| . 

We next discuss the calculation of /9-decay half-lives for Gamow- Teller decay and the re- 
lated problem of calculating /3-delayed neutron-emission probabilities. In our discussion of the 
model we use, unless otherwise stated, expressions and notation from the books by deShalit and 
Feshbach^^) and Preston and from our previous publications. ^'''^^) 
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2.5.1 /?- and /?+ decay 

The process of /3 decay occurs from an initial ground state or excited state in a mother nucleus 
to a final state in the daughter nucleus. For /3~ decay, the final configuration is a nucleus in 
some excited state or its ground state, an electron (with energy E^), and an antineutrino (with 
energy Ei,). The transition from the initial to the final state then involves an operator H, which 
is the weak-interaction Hamiltonian density. Once the operator H is known, the probability 
per unit time for emitting an electron with momentum between hke and ^(fee + dke) and an 
antineutrino with momentum between hk^, and h{ki, + dk^) is given by the well-known Golden 
Rule 

where Eq is the energy released in the decay. 

In the above expression one should sum over the spins of the final states and average over 
the initial spins. Our interest here is mainly to obtain the probability of decay to a specific final 
nuclear state /. To obtain this probability one must go through several lengthy steps. These 
steps are usually glossed over in discussions of these models, but one fairly extensive account 
of these steps is given in the book by Preston. ^'*) The final expression obtained through these 
steps for the total probability for decay to one nuclear state is 

^fi = '^^,Wfi?f^Z,R,e,) (8) 

where eo = Eo/mQC^, with rriQ the electron mass. For consistency with standard treatments of 
P decay we here use SI units. Moreover, |M^jp is the nuclear matrix element, which is also the 
/3-strength function. The dimensionless constant F is defined by 



where g is the Gamow- Teller coupling constant. There is a misprint concerning this quantity 
in the book^^) by deShalit and Feshbach, where in their Chapter 9, Eq. (2.11) the exponent 
is erroneously given as 2 instead of the correct value 3. The quantity f{Z, R, eo) has been 
extensively discussed and tabulated elsewhere. ^'^^^■'^) 

For the special case in which the two-neutron separation energy 6*211 in the daughter nucleus 
is greater than the energy Qp released in the decay, the probability for /3-delayed one-neutron 
emission, in percent, is given by 

Pin = 100 (10) 

0<Ef<Qf3 

where Ef is the excitation energy in the daughter nucleus and Si^i is the one- neutron separation 
energy in the daughter nucleus. We assume that decays to energies above always lead to 
delayed neutron emission. In the more general case where multiple-neutron emission is energeti- 
cally possible, but under the restriction that S'j/n in the daughter nucleus monotonically increase 
with increasing i/, we define 

PCk = 100 '^•"S"^^ (11) 

0<Ef<Qf3 
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where P^n probability of emitting v or more neutrons. The occupation probabihties Pa-u 

introduced above are then given by 

Pa = 100 -Pi> 

Pa-u = P.'n-Ptu+Dn (12) 

For some very neutron-rich nuclei S^n no longer monotonically increase with increasing zv, in 
which case the above formalism cannot be used. In the Table we present calculated values only 
when the above formalism is valid and otherwise put ". . . " in the columns for Pa-u- 

To obtain the half-life with respect to /3-decay one sums up the decay rates wji to the 
individual nuclear states in the allowed energy window. The half-life is then related to the total 
decay rate by 

In 2 

= (13) 

0<Ei<Qfi 

The above equation may be rewritten as 

„ n 27r3ln2 1 B 

(14) 



moc2 r2 ^ |M;,|2/(^,i2,eo) E \Mf^\^f{Z,R,eo) 

0<Ef<Qi3 0<Ef<Qp 

with 

moc^ i 

For the value of B corresponding to Gamow- Teller decay we use 

B = 4131 s (16) 

The energy released in electron emission is 

Qp- = = [M{Z, N) - M{Z + 1, iV - 1)] (17) 

whereas the energy released in positron emission is 

Qp+ = E^^ = [M{Z, N) - M{Z -1,N + 1)- 2mo] (18) 

The above formulas apply to f3~ and /J"*" decay. However, for calculating half-lives electron 
capture (EC) must also be considered. 



2.5.2 Electron capture 

The energy released in electron capture is 

Qec = eI;^ = [M{Z, N) - M{Z -l,N + 1)] - electron binding energy (19) 

so that 

E^''^ = Eq + 2mQ(? — electron binding energy (20) 

This shows that for some decays electron capture is possible whereas /J"*" decay is energetically 
forbidden. The total probability for decay to one nuclear state is again given by Eq. (|8|), where 
the final state / now refers to electron capture over all electron shells or to decay. 
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The total half-life with respect to /3+ and EC decay is given by 

In 2 

Tp = ^ (21) 

As pointed out above, the energies involved in the two terms in the sum differ by 2moc^ minus 
the electron binding energy and for some nuclear final states Wj^ may be zero (energetically 
forbidden) while w^f^ is not. 

To obtain an initial feel for these models we first studied the approximate relativistic expres- 
sions given by Preston. ^^) The results obtained by use of these expressions for f3~ decay are 
typically within 20% of those obtained in the more exact treatment by Gove and Martin, ^^) who 
have made extensive tabulations of /(Z, i2, eo). Despite the small differences we have obtained 
the computer code used to generate the tables of Gove and Martin and incorporated it into our 
programs. The results presented here have been obtained with this more accurate treatment. 

In the Table we present five quantities related to (3 decay, namely Pa, Pa-i, Pa-2i Qi3, 
and T^. The precise meaning of these quantities is as follows. When both /3"^ and (3~ decay 
are possible, we tabulate it in the columns for and T^. When neither /J"*" nor /3~ decay is 
possible, we tabulate ". . ." in the column for Qp, "/3-stable" in the column for Tp, and blank 
fields in the columns for Pa, Pa-i, and Pa-2- When only EC or f3~^ decay is possible, we 
tabulate Qec in the column for Q^, the calculated half-life with respect to Gamow- Teller decay 
for combined EC and f3~^ decay in the column for T^, and blank fields in the columns for Pa, 
Pa-1, and Pa-2- Finally, when only (3~ decay is possible, we tabulate Q^- in the column for 
Ql3, the calculated half- life with respect to Gamow- Teller f3~ decay in the column for T^, and 
the calculated occupation probabilities after /3-delayed neutron emission in the columns for Pa, 
Pa-1, and Pa-2- The electron binding energy has been neglected in the determination of Qec- 

To obtain more accurate values of and Pi,n, we have calculated the Q values that enter 
Eqs. ([Tl|), (13), and ( ^l] ) from experimental mass differences when all experimental masses that 



are required for the calculation are available and otherwise from calculated mass differences. The 
neutron separation energies that enter in the calculations of Pj^n in Eq. ( [ll|) are also obtained 
from experimental mass differences when available and otherwise from calculated mass differ- 
ences. Calculated deformations are always used, even when experimental data are available. 
However, in the astrophysical applications presented below, further use is made of experimental 
information, as is discussed in Sec. ||. 

2.6 Proton separation energies 

The one- and two-proton separation energies Sip{Z,N) and S2p{Z,N) are obtained from mass 
excesses through the differences 

SipiZ,N) = M{Z -l,N) + Mh- M{Z,N) 

S2p{Z,N) = M{Z -2,N) + 2Mn- M{Z,N) (22) 

2.7 a-decay properties 

The five heaviest known elements 107, 108, logMt, 110, and 111 were all identified from their 
a-decay chains, ^^~^'^) which limited their stability. The neutron number of the first identified 
isotope of 108 and Mt was A'^ = 157. Because a-decay chains provide very clear signatures of 
the nuclear species in the beginning of the decay chain and fission does not, it is likely that 
additional new nuclei discovered in the heaviest region will often be identified through their 
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a-decay properties. Models of a-decay properties are therefore highly useful for designing and 
interpreting experiments that explore the limits of stability of the heaviest elements. Obviously, 
one also needs to consider whether spontaneous-fission half-lives are significantly shorter than 
the a-decay half-lives. In that case spontaneous-fission would be the dominating decay mode 
and a decay might not be detected. 

The single most important quantity determining the a-decay half-life is the energy release 
Qa. In the heavy-element region an uncertainty of 1 MeV in Qa corresponds to an uncertainty 
of 10^^ for Qa ~ 7 MeV and to an uncertainty of 10^^ for ~ 9 MeV. The energy release 
Qa is obtained from the mass excesses through the difference 

Qa{Z, N) = M{Z, N) - M{Z -2,N -2)- M(2, 2) (23) 

The a-decay half-lives Ta presented in the Table are estimated by use of the Viola-Seaborg 
relationship ^) 

log(T„/s) = {aZ + 6)(Q„/MeV)-^/2 + (cZ + d) (24) 

where Z is the proton number of the parent nucleus. Instead of using the original set of constants 
suggested by Viola and Seaborg we use the more recent values 

a= 1.66175, b= -8.5166 
c= - 0.20228, d= - 33.9069 

that were determined in an adjustment taking into account new data for even-even nuclei.^) 
The uncertainties in the calculated half-lives due to this semi-empirical approach are far smaller 
than uncertainties due to errors in the calculated energy release. 

3 TABULATED RESULTS 

Deformed single-particle models provide the starting point for the calculations of nuclear ground- 
state masses and deformations, which were extensively discussed in our previous paper. ^) Since 
nuclear wave functions are also provided by these models, one may also use these models to 
determine electromagnetic moments and transition rates, /3-strength functions, /3-decay half- 
lives, and /3-delayed neutron-emission probabilities. 

The results of our calculations of many such nuclear properties for astrophysical applications 
are presented in the Table. To provide an overview of these results, we present in Figs. p]-p^ 
color diagrams of the calculated pairing gaps, neutron separation energies, proton separation 
energies, a-decay energy release and half-life, /3-decay half-life, and some related quantities. 
Beyond Z ~ 120 and N ~ 190 the calculated potential-energy surfaces on which Figs. p]~|l^ are 
based are very flat and the barrier with respect to fission is almost zero. Because the ground 
state is identified as the deepest local pocket in this fiat surface, one obtains rapidly fluctuating 
deformations, separation energies, and energy releases between neighboring nuclei. These results 
are of no physical significance, since the spontaneous-fission half-life of these nuclei will be much 
too short to allow experimental observation. 

3.1 Odd-nucleon spin and parity 

The most important constants in the folded- Yukawa single-particle model are the diffuseness 
and spin-orbit constants, which were determined '^^) in 1974 in the rare-earth and actinide re- 
gions from comparisons between calculated and experimental single-particle level orderings. The 
global nuclear-mass study ^^) in 1981 introduced a set of constants valid for the entire nuclear 
chart in terms of an expression for the spin-orbit strength that is linear in A = N + Z, with 
the expression fully defined by the previously determined values in the actinide and rare-earth 
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regions. This procedure is somewhat subjective because it is not based on exact comparisons 
between all available experimental data and calculations. Instead, one typically proceeds by cal- 
culating single-particle level diagrams as functions of deformation for several sets of constants, 
comparing their structure to a few selected nuclei, and choosing the set that gives the best 
agreement. 

Because we now have available nuclear ground-state shapes from our calculations of ground- 
state masses, we are in a position to compare calculated and experimental ground-state spins and 
parities in a well-defined manner, as shown in Figs. 17 and The only ambiguity concerns the 



comparison for nuclei calculated to be weakly deformed. We have chosen to base the comparison 
on spherical assignments if |e2| < 0.15 in the calculations. With this choice we obtain agreement 
in 446 cases and disagreement in 267 cases, corresponding to 63% agreement. When the ground- 
state energy is not minimized with respect to 63 and eg, we obtained in a previous calculation, 
as expected, slightly less favorable results: agreement in 428 cases and disagreement in 285 
cases, corresponding to 60% agreement. These results are not very sensitive to changes in the 
choice concerning when to use spherical assignments. In the present study, if we always choose 
the spherical assignment when this choice yields agreement with experimental data, we obtain 
agreement in 482 cases and disagreement in 231 cases, corresponding to 68% agreement, so that 
the improvement in the agreement is only 5%. 

The disagreements between the calculated and experimental spins and parities usually arise 
because several deformed or spherical levels lie very close together, making accurate calculations 
difficult. For magic numbers there is an almost stunning agreement, which, taken together with 
our analysis of the disagreements in other regions, makes it unlikely that a significantly better 
global set of constants can be found. The existing disagreements probably arise from residual 
interactions outside the framework of the single-particle model. 

3.2 Pairing gaps 

We first emphasize the most important results of our previous pairing study: ^) 

1. The preferred form of the pairing gap given by Eq. (^) lowers the rms deviation by about 
20% relative to the rms deviation obtained with the standard choice cj^fA for the average 
pairing- gap A or effective- interaction pairing gap Ag-. 

2. The Lipkin-Nogami pairing model yields an rms deviation that is 14% lower than the rms 
deviation in the BCS approximation. 

3. One cannot deduce the optimum constants for a microscopic pairing model by simplified 
macroscopic calculations. 

4. It is necessary to distinguish between several pairing-gap concepts, notably the average 
pairing gap A, the effective- interaction pairing gap Ac used as input in microscopic cal- 
culations, and the microscopic pairing gap A obtained as a solution to the BCS or Lipkin- 
Nogami pairing equations. 

5. The effective-interaction pairing-gap Ac does not depend explicitly on the relative neutron 
excess /. 

Our first four color figures represent results of theoretical calculations of pairing quantities 
in the Lipkin-Nogami approximation. In this model a pairing gap A and number-fluctuation 
constant A2 are obtained as solutions of the pairing equations. It is the sum A + A2, which we 
denote by Alnj that should be compared to odd-even mass differences. This quantity is shown 
in Fig. |l| for protons. The areas enclosed in jagged black lines are regions where experimental 
pairing gaps can be extracted from experimental masses by use of fourth-order finite-difference 
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Table 1: Errors of pairing-gap calculations. The experimental data base contains 1654 masses.^) 



LN Mass model 



Nucleons N^^^ /^th o"th crth;/i=o rms ^Uth dth o-th;/.=o rms 

(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) 



Neutrons 756 0.0220 0.1610 0.1627 0.1712 -0.0427 0.1929 0.1972 0.2088 
Protons 648 0.0455 0.1529 0.1596 0.1670 -0.0096 0.1728 0.1730 0.1994 
Total 1404 0.0328 0.1577 0.1613 0.1692 -0.0274 0.1846 0.1865 0.2045 



expressions. Because magic numbers and the N = Z Wigner cusp give non-smooth contributions 
to the mass surface and corresponding contributions to the finite-difference expressions, it is 
not possible to extract experimental pairing gaps from mass differences in certain regions near 
magic numbers and near N = Z. '^'^^) The proton pairing gap AlNp on the whole decreases 
with increasing A and with increasing Z. There is a decrease in the individual contribution 
Ap, which is not shown or tabulated in this paper, at magic proton numbers. This decrease is 
compensated for by a strong increase in A2p, which results in a smooth appearance of AlNp at 
magic proton numbers. 

The calculated neutron pairing gap AlNh is shown in Fig. |2[ The regions for neutrons where 
experimental pairing gaps can be extracted from odd-even mass differences are slightly different 
from those for protons shown in Fig. |l]. The individual contribution A^ decreases considerably 
near magic neutron numbers. This decrease is compensated for by a strong increase in \2n, so 
that their sum behaves relatively smoothly at magic neutron numbers, as seen in Fig. The 
pairing gap AlNh on the whole decreases with increasing A and with increasing N . There is no 
collapse at magic neutron numbers, in contrast to results based on the BCS approximation. 

As discussed in our previous pairing study, ^) one may also determine pairing gaps directly 
from odd-even mass differences based on theoretical masses. This type of theoretical pairing gap 
we denote by Athmass- There are several strong reasons to expect that it is more appropriate to 
compare Athmass to Aexp than to compare Aln to Aexp • The odd-even mass differences pick up 
any non-smooth contributions to the nuclear-mass surface, in addition to those represented by 
the pairing gaps. Although we have excluded regions near magic numbers from consideration, 
there are other gaps in the level spectra that are expected to give non-smooth contributions in 
the odd-even mass differences. These include the = 56 spherical subshell near the proton gap 
at Z = 40 and the A^ = 152 deformed gap near the proton gap at Z = 100 in the actinide region. 
Shape transitions also give rise to non-smooth components to odd-even mass differences. Such 
additional contributions to Aexp are expected near A^ = 88 in the transition to the deformed 
rare-earth region, for example. Figures ^ and Q clearly exhibit discrepancies related to these 
gaps and shape transitions that cannot automatically be interpreted as an inherent deficiency 
of the microscopic pairing model. 

These non-smooth contributions should be equally present in both Aexp and Athmass i Pro- 
vided that the theoretical calculation accurately describes the nuclear mass surface. Although 
neither Aexp nor Athmass represent the true pairing gap because of these non-smooth contribu- 
tions, they should cancel out in the discrepancy Aexp — Athmass- In Table || we present a study 
of the two discrepancies Aexp — Aln and Aexp — Athmass- In contrast to our expectations, we 
find better agreement between Aexp and Aln than between Aexp and Athmass- The total second 
central moment is 16% larger in the latter comparison than in the former. The reason that we 
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obtain such a large error for Athmass is that the deviation of the calculated mass surface from the 
true mass surface is sufficiently large to cancel the advantages of using Athmass as a theoretical 
pairing gap instead of Aln- In the transition regions between spherical and deformed nuclei 
near magic numbers, where we had expected such advantages, the error in the calculated mass 
surface is especially large, as can be seen in the comparison of measured and calculated masses 
in Fig. |l9. 



3.3 Total binding energy 

The total binding energy is equivalent to the atomic mass excess, which has been discussed 
extensively in our mass paper. ^) It is listed in the Table to facilitate applications that require 
binding energy rather than mass excess. 



3.4 Neutron separation energies 

One- and two-neutron separation energies are shown versus N and Z in Figs. The r-process 
takes place primarily in the region 2.0 MeV < 5in < 3.0 MeV for odd A^. Near = 50 and 
N = 82 the region of known nuclei extends into the r-process region. 

The discrepancy between experimental one-neutron separation energies obtained from mass 
differences and calculated one-neutron separation energies is shown in Fig. The biggest errors 
occur near magic numbers. There is a general decrease of the error as A increases. 



3.5 /5-decay properties 

A detailed knowledge of the low-energy part of the /3-strength function is essential for the 
calculation of many nuclear-structure quantities of astrophysical interest, such as the probability 
of /^-delayed neutron and proton emission, the probability of /3-delayed fission, and half-lives with 
respect to /3 decay. Experimentally it has been known for some time that the low-energy part of 
the /3-strength function exhibits a pronounced structure, ^5,36^ where the strength is collected in 
a few well-localized peaks. For nuclei that are spherical in their ground state there are usually 
only a few peaks within the window; for deformed nuclei the strength is more spread out, 
but still exhibits significant structure. 

Theoretically, these main features of the /3-strength function can be understood on the 
basis of an extreme single-particle model. The peaks in the strength functions correspond 
to transitions between specific single-particle levels. In the spherical case the levels are highly 
degenerate and spaced far apart, which gives rise to very few but strong peaks in the experimental 
strength function. For deformed nuclei the degeneracy is removed, allowing for significantly more 
transitions. Thus, compared to the spherical case, there are now more peaks in the experimental 
strength function, but the strength of each peak is lower. 

Although an extreme single-particle model explains the origin of the structure in the /3- 
strength function and the characteristic difference between strength functions associated with 
deformed and spherical nuclei, a more detailed description of the strength function requires 
the consideration of the residual pairing and Gamow- Teller interactions discussed above. The 
inclusion of these terms in the potential reduces the calculated strength in the low-energy part of 
the strength function to about 10% of what is obtained in an extreme single-particle model. 
Because pairing leads to a diffuse Fermi surface and consequently to some occupation probability 
above and to partially unfilled levels below the Fermi surface, there are decay channels open in the 
more refined model that are blocked in the extreme single-particle model. For deformed nuclei 
one often finds considerable strength for transitions between Nilsson levels whose asymptotic 
quantum numbers do not allow for any transition probability according to the GT selection 
rules. This occurs because the conventional, asymptotic quantum-number label gives the main 
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component of the wave function corresponding to the level, but the transition strength is due to 
small admixtures of wave functions with other asymptotic quantum numbers. Since we perform 
a full diagonalization of the single-particle Hamiltonian we account for these admixtures in our 
model. 

The /3-strength function is a sensitive measure of the underlying single-particle structure. 
Therefore, when calculating a /3-strength function, two conditions should ideally be fulfilled. 
First, the ground-state shape of the system of interest must be known. Second, the single-particle 
spectrum calculated at this shape must agree reasonably well with the experimental situation, 
especially for the levels closest to the Fermi surface. The calculation of a /3-decay half-life also 
requires the energy released in the decay, or equivalently the ground-state masses of the mother 
and daughter nuclei. We now use the folded- Yukawa single-particle potential, which has been 
applied to the calculation of nuclear masses, shapes, and other ground-state quantities for nuclei 
throughout the periodic system, i'33,37^ where these requirements are met globally. We noted in 
Sec. 3.1 that for spherical nuclei there was excellent agreement between calculated and measured 
ground-state spins and parities. For deformed regions disagreements occurred somewhat more 
often. However, when a disagreement occurs, a level with the correct spin and parity is often 
calculated to be very near the last occupied orbital. If this situation occurs in the daughter 
nucleus there is a rather small effect on the calculated /3-strength function. If, however, in a 
decay from an odd-even or odd-odd nucleus the unpaired proton or neutron is in an incorrect 
level the effect can be more significant. Therefore, it is reasonable to expect characteristic 
differences in the discrepancies between calculated and experimental half-lives in odd-even and 
even-even decay. Surprisingly, as will be discussed below, we see no such differences. 

To illustrate some major features associated with /3-decay we show in Figs. 20 and ^ cal- 
culated /^-strength functions for the spherical nucleus ^^Rb and the deformed nucleus ^^Rb. 
Although we obtain a fairly large deformation for ^^Rb in our calculation of nuclear ground- 
state masses and deformations^) we have performed the calculation shown in Fig. 20 for a 
spherical shape, which is consistent with experimental information on this nucleus, so that we 
can study /3-strength-function features associated with a sudden onset of deformation as the 
neutron number increases from = 58 to = 62. As discussed in more detail in Ref. ^®), 
the single-particle and pairing properties are evaluated for the appropriate vacuum nucleus. In 
addition, the Qp value corresponds to the decay indicated and the Sy^ and Py^ are evaluated for 
the daughter of the decay. A characteristic difference between /^-strength functions of spherical 
and deformed nuclei is that the former contain only a few strong peaks in the low-energy region, 
whereas the latter contain many more peaks of smaller size. Therefore, the likelihood that one 
or several of them will occur close to zero energy is much larger in the deformed case than in 
the spherical case. Because the transition rate is proportional to Sp{Ex — Qp)^, where Ex is the 
excitation energy of the daughter state, there is therefore usually a characteristic drop in the 
/3-decay half-life at the transition from a spherical to a deformed system. 

It is not our aim here to make a detailed analysis of each individual nucleus, but instead 
to present an overview of the model performance in a calculation of a large number of (3- 
decay half-lives. In Figs. 22|-27 we compare measured and calculated /3-decay half-lives for 
nuclei throughout the periodic system. The experimental half-lives are from a compilation by 
Browne. We present results with respect to /3~ decay and also to decay and electron 
capture. To avoid lengthy constructions we will in our discussion in this section usually not 
distinguish between /J"*" decay and electron capture and somewhat inexactly take (5/3 to mean 
the maximum energy release in the decay. Because the calculated pairing gap affects o-o, o-e, 
and e-e decays differently we analyze these decays separately to differentiate between effects due 
to pairing and due to other causes. 

We have limited the comparison to nuclei whose experimental half-lives are shorter than 
1000 s. Because the relative error in the calculated half-lives is more sensitive to small shifts in 
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the positions of the calculated single-particle levels for decays with small energy releases, where 
long half-lives are expected, one can anticipate that half-life calculations are more reliable far 
from stability than close to /3-stable nuclei. To address the reliability in various regions of nuclei 
and versus distance from stability, we present in Figs. 22-|27 the ratio T/3,caic/7/3,exp versus the 
three different quantities iV, T^^exp; and Qp. 

The few cases that lie outside the scale of the figures have been taken into account in the 
error analysis presented below. These cases are associated with long experimental half-lives; 
there is only one case outside the scale of the figure with an experimental half-life less than 
10 s. In such cases first-forbidden transitions, which we do not take into account at this stage, 
can be expected to make a considerable contribution to the decay rate. Therefore, one should 
not conclude from these relatively few very large deviations that the calculated Gamow- Teller 
/3-strength function is in serious disagreement with experimental data. 

Before we make a quantitative analysis of the agreement between calculated and experimental 
half-lives we briefly discuss what conclusions can be drawn from a simple visual inspection of 
Figs. p^27. In Figs. 22 and 25 the quantities I/3,caic/r/3,exp are plotted as functions of neutron 



number N . There are no systematic trends with N . For (5 decay of even-even nuclei the 



calculated half-lives are somewhat too long on the average. In Figs. |23| and |26| the quantities 
7/3,caic/T/3,exp are plotted as functions of the experimental half-life T^^cxp- As a function of this 
quantity, one would expect the average error to increase as TJg^exp increases. This is indeed 
the case for j3~ decay, but not for /?+ decay, except for odd-odd decays — an unexpected result. 



In Figs. 24 and 27 the quantities J/3^caic/T'/3,exp are plotted as functions of Qp with the aim of 
showing how the average error increases as Q/3 decreases. It is obvious that errors in the location 
of the peaks in the calculated strength function have a larger effect on the calculated half-lives 
for small values of Q/j than for larger ones. Indeed, we see a fairly clear increase in the scatter 
of the points in Fig. ^ as Qp decreases. However, for decay shown in Fig. ^ the only 
correlation occurs for 0-0 decay. 



In a visual inspection of Figs. 22-27 one is left with the impression that the errors in our 
calculation are fairly large. However, this is partly a fallacy, since for small errors there are 
many more points than for large errors. This is not clearly seen in the figures, since for small 
errors many points are superimposed on one another. To obtain a more exact understanding of 
the error in the calculation we therefore perform a more detailed analysis. 

One often analyzes the error in a calculation by studying a root-mean-square deviation, 
which in this case would be 

1 " 2 

<7rms^ = — '^{Tp^e^p — T^,calc) (25) 
i=l 

However, such an error analysis is unsuitable here, for two reasons. First, the quantities studied 
vary by many orders of magnitude. In our case the variation is more than ten orders of mag- 
nitude, from the millisecond range to years and beyond. Second, the calculated and measured 
quantities may dijjerhy orders of magnitude. We therefore study the quantity log(T^^caic/7'/3,exp)) 
which is plotted in Figs. 22-27, instead of (TJg^exp — 7/3,caic)^- 

To facilitate the interpretation of the error plots we consider two hypothetical cases. As 
the first example, suppose that all the points were grouped on the line T/3^caic/T/3,exp = 10. 
It is immediately clear that an error of this type could be entirely removed by introducing a 
renormalization factor, which is a common practice in the calculation of /?-decay half-lives. We 
shall see below that in our model the half-lives corresponding to our calculated strength functions 
have about zero average deviation from the calculated half-lives, so no renormalization factor is 
necessary. 

In another extreme, suppose half the points were located on the line TJg,caic/?/3,cxp = 10 and 
the other half on the line 7/3,caic/7/3,exp = 0.1. In this case the average of log(T^^caic/7)3,cxp) 
would be zero. We are therefore led to the conclusion that there are two types of errors that are 
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of interest to study, namely the average position of the points in Figs. 22-21 ^ which is just the 



average of the quantity log(T/3^caic/7/3,oxp), and the spread of the points around this average. To 
analyze the error along these ideas, we introduce the quantities 

f = T/3^calc/T/3,exp 

n = logio(?-) 

1 



Mr 



10 



n 



a. 



a, 



1 " 



1=1 



1/2 



(26) 



where M^j is the average position of the points and is the spread around this average. The 
spread CTrj can be expected to be related to uncertainties in the positions of the levels in the 
underlying single-particle model. The use of a logarithm in the definition of n implies that these 
two quantities correspond directly to distances as seen by the eye in Figs. |2^-p7|, in units where 
one order of magnitude is 1. After the error analysis has been carried out we want to discuss its 
result in terms like "on the average the calculated half-lives are 'a factor of two' too long." To 
be able to do this we must convert back from the logarithmic scale. Thus, we realize that the 
quantities M^^ and 0"^^'^ are conversions back to "factor of" units of the quantities M^^ and fT^j , 
which are expressed in distance or logarithmic units. 

In Table ^ we show the results of an evaluation of the quantities in Eq. ( [2^ ) for f3~ decay 
of nuclei with ?/3,exp ^ 1000 s. For long half-lives one expects forbidden decay to dominate. 



Table 2: Analysis of the discrepancy between calculated and measured /3 -decay half-lives shown 



in Fig. 23 
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Mr, 








T^max 
/3,exp 


o-o 


29 


-0.23 





59 


0.46 


2.91 
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o-e 


35 


-0.23 





59 


0.42 


2.64 
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e-e 


10 


0.58 


3 


84 


0.49 


3.08 
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o-o 


59 


-0.12 





76 


0.95 


8.83 


10 
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85 


-0.11 





78 


0.68 


4.81 


10 


e-e 


34 


0.40 
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50 


0.62 


4.13 


10 


o-o 


88 


0.37 


2 


33 


1.69 


49.19 


100 


o-e 


133 


0.05 


1 


11 


0.98 


9.45 


100 


e-e 


54 


0.42 


2 


61 


0.68 


4.75 


100 


o-o 


115 


0.54 


3 


50 


1.86 


72.02 


1000 


o-e 


194 


0.44 


2 


77 


1.85 


71.50 


1000 


e-e 


71 


0.84 


6 


86 


1.77 


58.48 


1000 
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We find that in the o-o and o-e cases the value of M^^ increases from values somewhat below 
1 to values somewhat above 1 as the maximum half-life, T^f^„, considered increases from 1 s 
to 1000 s. When cases corresponding to long half-lives contribute to the error estimates, first- 
forbidden decays would be expected to increase the /3-decay rates. In our model we neglect, 
at this stage, first-forbidden decay. Therefore, one can expect that as T™^^ increases, our 
calculated half-lives will become increasingly longer than experiment. This is indeed the trend 
in Table ^. Therefore, it is preferable to use only nuclei with short half-lives to determine 
if any renormalization of the /3-strength function is needed to reproduce on the average the 
experimental half-lives. However, for a short half-life cutoff, for example 1 s, there are not enough 
data points for a reliable determination. Slightly higher half-life cutoffs indicate M}^ ~ 1, except 
for the e-e case. However, there are only a few data points for the e-e case, so the higher M}^ 
obtained in this case may be unreliable, especially in view of the results obtained below for (5'^ 
decay. 

In Table |3| we show the results of an evaluation of the quantities in Eq. ( ]26| ) for decay and 
electron capture of nuclei with 7/3,exp < 1000 s. Above we argued in the case of j3~ decay that 
as nuclei with longer half-lives are included, one can expect the calculated half-lives to be too 
long because the rates due to first-forbidden decay are not included. However, we see instead 
a slightly opposite trend here. But from arguments similar to those used in the discussion of 
f3~ decay we conclude that in /J"*" decay our results are consistent with zero average deviation 
between calculated and experimental half-lives. Thus, the error in the calculation is given by 
(T^l^ only. The behavior of al^ as TJ^^^ increases is similar for both j3~ and /?"^ decay. Just as 
when the average deviation is determined, one should not consider decays with too-long half- 
lives when determining a^^ . On the other hand, it would be desirable to base the determination 
of a}!^ on a large set of data. Both of these conditions cannot be well-fulfilled simultaneously. 
However, if we base our estimate of the model error cj^J^ by weighting most heavily half- lives in 
the shorter range, we find that the model error is about a factor of four. 

Certain additional assumptions and approximations have been made to extend the odd-^ 

Table 3: Analysis of the discrepancy between calculated and measured /3"'"-decay and electron- 
capture half-lives shown in Fig. |2^. 
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QRPA formalism to the odd-odd case. If these simphfications were inadequate, one would expect 
to see a larger 0"^ in the o-o case than in the o-e case. No clear effect of this type is seen in 
Table ^ for the shorter half-lives but does develop for very long experimental half-lives. 

We have already mentioned that particle-particle correlations are not taken into account in 
our QRPA treatment. These correlations are expected to strongly suppress /3+ transitions, 2i>39) 
and as a consequence lead to longer half-lives than those obtained in a model that does not 



consider these correlations. However, our results in Figs. 22-27 and Tables g and y show that 
in a standard treatment with xgt = 23 MeV/j4 and no renormalization coefficient, we obtain 
/3+-decay rates that actually agree better with experiment than do our calculated /3~-decay 
rates. 

3.6 Proton separation energies 

One- and two-proton separation energies are shown versus and Z in Figs. |9-11. For odd Z the 
region of known nuclei extends into the region where proton emission is energetically allowed in 
several places, but especially just above N = 82. However, the Coulomb barrier severely inhibits 
proton emission from the ground state when the energy released is small. Proton emission is 
therefore more readily observed as delayed emission following (3 decay. '^'^'^^ ) 

The discrepancy between experimental one-proton separation energies obtained from mass 



differences and calculated one-proton separation energies is shown in Fig. 12. The biggest errors 



occur near the magic numbers. There is a general decrease of the error as A increases. 
3.7 ci-decay properties 

The calculated energy release and associated half-life with respect to a decay are plotted in 
Figs. ^ and 14, respectively. 

In 1989 Miinzenberg et al.'^'^) compared the a-decay energy release along the N = 154 
and 155 isotonic lines to predictions of the FRLDM (1988).^^) In Fig. |2^ we make a similar 
comparison of measured data to predictions of the current FRDM (1992).^) These results 
based on the current FRDM (1992) agree much better with the measured values than do the 
corresponding results calculated from our older mass model in 1989 by Miinzenberg et al.,'^'^) 
and also much better than the results obtained with the 1991 version of the FRDM. ■*^) The 
improvement is due partly to the inclusion of the eg shape degree of freedom. 

Prom the Table we find that the nucleus ■^''^110 has a calculated a-decay half-life of 71 ms. 
The nuclei -^^^110 and ■^^'^110 in the center of the superheavy island have calculated a-decay half- 
lives of 4 y and 1565 y, respectively. With the 1991 version of the FRDM we obtained calculated 
a-decay half-lives of 16 ms, 161 y, and 438 y for '^'^'^110, ^^^110, and ^^^110, respectively. Our 
current results are different from these older ones for two reasons. First, we now use the new 
mass model FRDM (1992). Second, we now use a new set of constants in Eq. (p4[). For ^^^110 
the value of Qa is 7.36 MeV in our current model, but only 6.95 MeV in the FRDM (1991). 
For ^^^110 the difference in Qa values is only 0.04 MeV, and the difference in the calculated 



half- lives is in this case due almost entirely to the new set of constants used in Eq. (24). 

In the introduction to Sec. ^ we pointed out that because the calculated potential-energy 
surfaces for some very proton-rich nuclei are very flat and contain multiple minima, the ground- 
state deformation may exhibit large fluctuations between neighboring nuclei. For example, we 
see in Ref. ^) that the calculated quadrupole deformation of ^^'^128 is €2 = 0.0 whereas €2 = 0.392 
is obtained for '^^^128. Such deformation changes are accompanied by non-smooth changes in 
calculated energy releases and associated half-lives. An a-decay half-life of lO"*^'^^ s is obtained 
for ^^^128, whereas the calculated a-decay half-life for '^^^128 is greater than lO^'' s. Because 
of very short spontaneous-fission half-lives for such proton-rich nonspherical nuclei, we do not 
expect these nuclei to be observable despite their calculated long a-decay half- lives. 
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4 EXTRAPOLATEABILITY 

Theories are often characterized by their rms error with respect to experimental data points. 
However, this is often unsuitable because the rms error contains contributions from errors in the 
experimental data points and is therefore always larger than the intrinsic error of the theory. 
On the other hand, it is possible to determine the intrinsic error of a theory in a way that 
contains no contributions from experimental errors. It is most natural to characterize the 
deviation between measured data points and theoretical calculations by two quantities: (1) the 
mean deviation /ith of the theory from the data points and (2) the standard deviation cJth of 
the theoretical results about this mean. Alternatively, one may characterize the error of 

a theory by a single number, the second central moment of the error term of the theory, which 
we denote by o"th;/i=o- This measure is similar to the rms deviation, but without the deficiency 
of contributions from experimental errors. We use these four error measures in our discussions 
of model accuracies below. 

After we completed our mass model, ^) a new compilation of experimental ground-state 
masses ^'^) has become available. We are therefore now able to compare our FRDM (1992) 
predictions to data that were not taken into account when the model constants were determined, 
instead of studying the model reliability through simulation. In Table |^ we give errors for the 
FRDM(1992) for ground-state masses and six separation energies and energy releases for both 
the region in which the model constants were determined and in the new region of nuclei. 

We see in Table |^ that the error o"th;/^=o foi' iisw masses, 0.642 MeV, is 4% smaller than the 
error 0.669 MeV in the region where the model constants were determined. Thus, we conclude, 
as was done in our mass paper ^) through simulations, that the model is very well-behaved as 
one moves away from stability to new regions of nuclei. Indeed, at this stage there is no sign 
that the model diverges. It is also very clear that it is not possible to use the rms deviation as 
a proper measure of model error. It is 7% larger in the new region of nuclei, although we have 
just shown that the true model error is actually 4% smaller in the new region of nuclei than in 
the region where the model constants were determined. This is because the experimental errors 
in the new region of nuclei are quite large and contribute significantly to the rms error, whereas 
they are normally quite small in the region of previously known nuclei. 

In the absence of correlations in the model error between neighboring masses, the model error 
for the separation energies and energy releases would be \/2 times the mass error. However, 



Table 4: Theoretical errors and extrapolateability of the FRDM (1992). 





Orig 


inal nuclei 


(1989) 




New nuclei {191 


B9-1993) 






^nuc 




rms 


Vnuc 




CTth 


CTth;/^=0 


rms 






(McV) 


(MeV) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 




1654 


0.669 


0.681 


217 


0.074 


0.638 


0.642 


0.730 


<S'ln 


1464 


0.393 


0.411 


195 


0.033 


0.525 


0.526 


0.624 


<S'2n 


1403 


0.539 


0.555 


198 


0.026 


0.575 


0.576 


0.684 


Qp- 


1353 


0.488 


0.507 


216 


0.006 


0.559 


0.559 


0.647 


Sip 


1400 


0.381 


0.397 


210 


0.068 


0.576 


0.581 


0.672 


S2p 


1314 


0.493 


0.509 


193 


0.026 


0.585 


0.586 


0.663 


Qa 


1450 


0.625 


0.641 


224 


0.021 


0.675 


0.675 


0.772 
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Table 5: Comparison of theoretical errors and extrapolateability of the FRDM (1992) with those 
of the ETFSI-1 (1992) model. 



Original nuclei (1989) New nuclei (1989-1993) 







FRDM 


ETFSI-1 






FRDM 




ETFSI-1 






^ * nuc 




rms 




rms 


*nuc 


Mth 




rms 




O'th;/i=0 


rms 






(MeV) 


(MeV) 


(MeV) 


(MeV) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


Mth 


1540 


0.607 


0.615 


0.733 


0.742 


213 


0.070 


0.645 


0.716 


0.015 


0.809 


0.866 


'S'ln 


1356 


0.335 


0.346 


0.514 


0.521 


190 


0.033 


0.525 


0.609 


0.044 


0.612 


0.681 


'S'2n 


1301 


0.477 


0.482 


0.501 


0.508 


192 


0.028 


0.580 


0.659 


-0.052 


0.576 


0.659 


Qp- 


1259 


0.425 


0.436 


0.679 


0.684 


212 


0.001 


0.560 


0.646 


-0.003 


0.741 


0.791 




1296 


0.324 


0.339 


0.466 


0.472 


205 


0.073 


0.581 


0.660 


0.113 


0.608 


0.709 


5'2p 


1217 


0.434 


0.444 


0.520 


0.531 


189 


0.030 


0.589 


0.662 


0.082 


0.586 


0.674 




1342 


0.541 


0.546 


0.520 


0.529 


217 


0.018 


0.670 


0.762 


0.004 


0.618 


0.699 



Table ^ shows that the separation-energy and energy-release errors are instead somewhat smaller 
than the mass-model error. This is due to the correlation in the mass-model error for neighboring 
nuclei. The separation-energy and energy-release errors increase in the new region of nuclei, 
whereas the mass- model error does not. This indicates that the correlation between mass- model 
errors for neighboring nuclei decreases in the new region. Since a mass model should not be 
judged on the amount of error correlations between nearby nuclei, one cannot gain significant 
insight about a model by comparing the error of the separation energies and energy releases in 
the region where the model constants were determined to the errors in a new region of nuclei. 
However, because many of these quantities are used in astrophysical calculations, there is a 
practical need to know the errors of these quantities, and it is for this purpose that we provide 
them. 

In Table | we compare the errors of the FRDM (1992) and the 1992 extended Thomas- 
Fermi Strutinsky-integral model (version 1) of Aboussir, Pearson, Dutta, and Tondeur ^'''^^) in 
the region of nuclei where the ETFSI-1 (1992) model constants were determined and in the 
region of new nuclei that are included in the ETFSI-1 (1992) calculation. The region of nuclei 
considered in the ETFSI-1 (1992) model is slightly different from the region considered in the 
FRDM (1992). Its constants were also determined from an adjustment to a slightly earlier mass 
evaluation. ) In Table || the quantities pertaining to the FRDM (1992) are also evaluated for 
this more limited region. Therefore, the results for the FRDM (1992) in Table | are slightly 
different from those in Table ^ The finalized ETFSI-1 (1992) mass table, as it appears in 
Ref. ^^), differs slightly from an earlier table circulated privately before publication. We use here 
the finalized, published version. 

The comparison in Table | between the extrapolateability of the FRDM (1992) and the 
ETFSI-1 (1992) model shows that for these particular regions the FRDM (1992) error increases 
by 6% and the ETFSI-1 (1992) model error increases by 10% in the new region relative to the 
region where the ETFSI-1 (1992) model constants were determined. Thus, we conclude from 
Tables | and I that the FRDM (1992) is somewhat more stable as one moves away from stability 
than is the ETFSI-1 (1992) model. 

In Table | we compare the extrapolateability between the FRDM (1992/1977) and the 1976 
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Table 6: Comparison of theoretical errors and extrapolateability of the FRDM (1992/1977) with 
those of the mass formula of von Groote et al. (1976). 



Original nuclei (1977) New nuclei (1977-1989) New nuclei (1989-1993) 



FRDM V. Groote FRDM v. Groote FRDM v. Groote 





Nnuc 


O'th;/x=0 


Cth;fi=0 


iVnuc 


Mth 


O'th;fj.=0 




O'th;fj.=0 


iVnuc 


^J■th 


O'th;/i=0 




O'th;p=0 






(MeV) 


(MeV) 




(MeV) (MeV) 


(MeV) (MeV) 




(MeV) (MeV) 


(MeV) (MeV) 


Mth 


1323 


0.671 


0.629 


351 


0.004 


0.686 


0.612 


1.154 


217 


0.030 


0.668 


0.808 


1.284 


'S'ln 


1139 


0.391 


0.383 


335 


-0.032 


0.406 


-0.107 


0.431 


195 


0.024 


0.525 


-0.095 


0.541 


5*211 


1092 


0.508 


0.480 


320 


-0.053 


0.632 


-0.171 


0.650 


198 


0.007 


0.576 


-0.213 


0.609 


Qp- 


1052 


0.467 


0.516 


314 


0.077 


0.552 


0.193 


0.737 


216 


0.031 


0.562 


0.178 


0.633 


Sip 


1105 


0.361 


0.412 


305 


0.039 


0.463 


0.102 


0.548 


210 


0.083 


0.585 


0.157 


0.619 


S2p 


994 


0.478 


0.501 


331 


0.049 


0.546 


0.160 


0.691 


193 


0.056 


0.591 


0.217 


0.738 




1144 


0.627 


0.525 


321 


-0.005 


0.630 


-0.001 


0.669 


224 


0.011 


0.675 


0.024 


0.674 



mass formula of von Groote, Hilf, and Takahashi. ^'') In this comparison we have adjusted the 
constants of the FRDM to similar masses ^^) as those considered by von Groote et al. in the 
determination of their constants. Clearly, the mass formula of von Groote et al. diverges severely 
in new regions of nuclei. 

In Figs. 29-37 we plot the deviations between calculated and experimental masses, one- 
neutron separation energies, and energy releases for /3~ decay for these three models. Tables ^|-^ 
and Figs. 35-37 are based on all Qfj- values, both positive and negative. Neither the FRDM 
(1992) nor the ETFSI-1 (1992) model shows large divergences in the new region, but strong 
systematic deviations occur in the model of von Groote et al. 

We have also made a limited study of the extrapolateability of the recent 1994 Thomas-Fermi 
model of Myers and Swiatecki. ^^'^^) In this model, the macroscopic energy is calculated for a 
generalized Seyler-Blanchard nucleon-nucleon interaction by use of the original Thomas-Fermi 
approximation. For A^, Z > 30 the shell and pairing corrections were taken from the 1992 finite- 
range droplet model, and for A^, Z < 29 a semi-empirical expression was used. The constants of 
the model were determined by an adjustment to the ground-state masses of the same 1654 nuclei 
with N,Z > 8 ranging from ^^O to ^^^106 whose masses were known experimentally in 1989 
that were used in the 1992 finite-range droplet model. The theoretical error corresponding to 
these 1654 nuclei is 0.640 MeV. The reduced theoretical error relative to that in the 1992 finite- 
range droplet model arises primarily from the use of semi-empirical microscopic corrections in 
the region N, Z < 29 rather than microscopic corrections calculated more fundamentally. The 
theoretical error for 217 newly measured masses is 0.737 MeV, corresponding to an increase 
of 15%. This model therefore extrapolates to new regions somewhat less well than does either 
the FRDM (1992) or the ETFSI-1 (1992) model. The deviations between experimental and 
calculated masses for these 217 new nuclei are shown in Fig. 

In the Wapstra-Audi mass evaluations ^~^'*^^ ) there usually are listed not only measured 
masses, but also masses estimated from systematic trends. It is unfortunate that some theoretical 
mass studies have not observed the difference between measured masses and masses given by 
systematic trends, and included also the latter in the data set to which the model parameters 
were adjusted. The masses given by systematic trends have to be considered on the same basis 
as masses given by other models. We may therefore study the reliability of the Wapstra-Audi 
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systematic-trend model by use of the methods apphed above to other models. We consider the 
masses given by systematic trends in the 1989 midstream evaluation^) as model masses and 
compare them to the new masses determined in the 1993 mass evaluation. ^''^) The result is 
shown in Fig. 39. The systematic-trend model error in the new region is 0.404 MeV, which 
is smaller than the error of the FRDM (1992). However, the systematic-trend model does not 
provide masses very far from known nuclei, and indeed many of the newly measured masses 
were not predicted by the 1989 systematic-trend model. Also, in this case it is not possible to 
provide a ratio between the error for nuclei in the new region and the error in the known region. 

The discovery in late 1994 of the new element Z = 111^") allows us to test model extrap- 
olateabilities in the region of large proton numbers. In Fig. Eo| we compare the experimental 
results for the a-decay chain of the heaviest known element 111 with predictions obtained 
in the FRDM (1992), the ETFSI-1 (1992) model ^^-"^S), and the 1992 fermion dynamical sym- 
metry model of Han, Wu, Feng, and Guidry. ^■*) Clearly the predictions of the FRDM (1992) 
agree much better with the new experimental data than do those of the other two models. The 
FRDM (1992) error for this chain is in fact considerably smaller than the error in the region 
where the model constants were adjusted, whereas the errors for the other two models, especially 
the FDSM (1992), are considerably larger than in the regions where the constants of these two 
models were adjusted. Because the FRDM (1992) Qa error is about 0.5 MeV for heavy nuclei in 
the known region, one should not consistently expect the exceptionally good agreement present 
in Fig. 40 for all new a-decay chains that are discovered in the heavy region. However, the 
agreement between experiment and predictions of the FRDM (1992) seen in this figure, plus 
better-than-expected agreement between experiment and predictions of the FRDM (1992) and 
the decay chains of another recently discovered element with proton number Z = 110, ^^) confirm 
the conclusion reached earlier by simulation^) that the FRDM (1992) is reliable to its stated 
accuracy as we move away from the region of known elements towards the superheavy region. 



5 ASTROPHYSICAL APPLICATIONS 

Nuclear physics and astrophysics share several common themes: (1) the nuclear reactions that 
are responsible for nucleosynthesis and isotopic abundance patterns in nature, (2) the energy 
sources of stellar events via static or explosive nuclear burning, and (3) the behavior of nuclear 
matter at and beyond nuclear densities, including the equation of state for type-H supernova 
(SN H) explosions and neutron stars. 

Here, we focus on explosive conditions in which large numbers of free protons or neutrons 
are available. They can lead to the synthesis of nuclei far from stability either via rapid proton 
capture and f3~^ decay (rp-process) in novae and X-ray bursts or rapid neutron capture and 
P~ decay (r-process) in SN H. In both cases, theoretical studies require reaction rates and 
other nuclear properties of unstable nuclei which, to a large extent, are not known experimen- 
tally. Hence, a general understanding of their nuclear-structure properties can be obtained only 
through theoretical means. 

Because several different nuclear quantities are needed in rp- and r-process calculations, in 
the past it was not possible to obtain them all from one source. Taking them from different 
sources, however, raises questions of consistency. In such mixed-model calculations, ^^) although 
often performed due to the lack of a unified approach, occasionally nuclear structure signatures 
may vanish or artificial effects may occur, thus strongly limiting their predictive power far from 
stability. Therefore, attempts to use this approach to identify the sites of rp- and r-processes 
may lead to erroneous conclusions. Consequently, it is of great importance to provide nuclear- 
structure properties based on a single, unified theoretical framework within which all quantities 
of interest can be obtained. We illustrate this principle below with an example. 
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5.1 The rp-process 

In a recent paper van Wormer et al.^^) have evaluated in detail the nuclear-reaction sequences 
for conditions in explosive hydrogen burning at temperatures beyond 10^ K, corresponding to 
the rp-process, a sequence of rapid proton-capture reactions and decays passing through 
proton-rich nuclei. Such processes typically occur when hydrogen fuel is ignited under highly 
degenerate conditions in explosive events on the surface of compact objects like white dwarfs 
(novae) and neutron stars (X-ray bursts). Hydrogen burning at high temperatures such as 
2 X 10^ K may also occur in a stable fashion in so-called Thorne-Zytkow objects. Such objects 
are expected to result from the merging of a neutron star and a main-sequence star in a binary 
system. An rp-process occurs then at the base of a fully convective envelope just above the 
neutron core, giving the object the appearance of a red supergiant. 

As pointed out by some authors, ^^'^'') the rp-process cannot be explained in a simple and 
clean fashion in terms of a (p,7) ^ (TjP) equilibrium in all isotopic chains. Under such circum- 
stances, the nuclear data needed for a theoretical description would be limited to nuclear masses 
and corresponding proton separation energies and /3"'"-decay half-lives. However, the Coulomb 
barriers in charged-particle captures lead to a cycle pattern with capture and decay time scales 



of similar size. As outlined in Fig. |4lj, the rp-process at low temperatures is dominated by 
two successive proton captures, starting out from an e-e nucleus, a 13'^ decay, a further proton 
capture into an even-Z nucleus, another /J"*" decay, and a final (p,a) reaction close to stability. 

The progress of the rp-process towards heavier nuclei depends on the leakage ratio (p,7) / (p,a) 
into the next cycle. Increasing temperature makes it possible to overcome Coulomb barriers 
and extend the cycles to more proton-rich nuclei, which permits additional leakage via proton 
captures competing with long (5 decays. Beyond 3 x 10® K, all cycles break open and a complete 
rp-pattern of proton captures and /? decays is established, which may reach beyond ^^Ni. ^6,57^ 

As a consequence, along with the separation energies and half-lives that are required for 
an equilibrium assumption, proton and a-capture rates for unstable nuclei are also required as 
nuclear-physics input into the reaction network. In current calculations the majority of these 
reaction rates are not based on measured cross sections, but are usually approximated by simple 
statistical- model calculations.^®) Even such statistical approaches require a number of input 
quantities, such as (1) ground-state spins and parities of target, compound, and final nuclei, (2) 
reaction energy releases, (3) realistic nucleon-nucleus and a-nucleus optical potentials, (4) giant- 
dipole resonance energies and widths, and (5) level-density models. It is desirable to improve 
the accuracy of current level-density models by utilizing shell and pairing corrections of modern 
mass models. 

5.2 The r-process 

The r-process was inferred ^^'^'^) from the observation of characteristic peaks in the abundance 
curve of /3-stable nuclei. The peak maxima are located 4, 4, and 9 neutrons below the magic 



neutron numbers 50, 82, and 126, respectively. In Fig. 16 some basic features of the r-process 
are outlined. The position of the r-process line depends on nuclear-structure properties and the 
stellar conditions under which it occurs, in particular the temperature, density, and duration 
of the neutron flux. The line plotted in the figure corresponds to a neutron separation energy 
of 2.4 MeV. The r-process is dynamic and may be located anywhere in the region 1.5 MeV 
< Sin < 3.0 MeV, depending on the stellar conditions. Magic neutron numbers play a special 



role in the r-process as is partially seen in Fig. 16, where the r-process line has a kink at each 



magic neutron number. In addition, nuclei pile up at these magic neutron numbers because 
of the long /3-decay half-lives and sudden lowering of the neutron-capture cross section that 
occurs at magic neutron numbers. After freezeout of the neutron flux, these nuclei {3 decay back 
towards the line of /3 stability. It is easy to see from the figure that the peaks in the abundance 
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curve are related to a concentration of neutron-magic nuclei far from (3 stability. In fact, this 
effect has been used to put rough constraints on nuclear masses far from /? stability and to rule 
out nuclear mass models with predictions far from these constraints. The decay back to the line 
of P stability is also influenced by /3-delayed neutron emission. 

Since the pioneering work of Burbridge et al. and Cameron ^'^) the rapid neutron capture 
process has been associated with explosive environments with high temperature T > 10^ K and 
high neutron-number density cm Under such conditions, the neutron-capture time 

scales of heavy nuclei are so short that within about 10~^ s highly neutron-rich nuclei can be 
produced up to 15 to 30 mass units away from (3 stability with neutron separation energies 
^in ~ 1.5 to 3.0 MeV. Neutron captures are not hindered by increasing Coulomb barriers, 
in contrast to the case for the rp-process in the previous section. Magic neutron numbers are 
encountered for smaller mass numbers A than in the valley of stability, which, after freezeout 
and (3 decay, shifts the observed solar-system r-process abundance peaks below the s-process 
peaks. However, besides this basic understanding, the history of r-process research has been 
quite diverse in suggested astrophysical scenarios, as well as with the required size of the nuclear 
network. 

The observed isotopic r-abundances A^r,© shown in Fig. 16 are the result of successive neu- 
tron captures along the r-process path and /?~ decay back to stability, thus depending, apart 
from stellar parameters, on a variety of nuclear properties of nuclei with extreme N/Z ratios. 
In the general case, among nuclear-physics quantities, ground-state masses and corresponding 
Qp values, neutron-separation energies S'ln, /3~-decay half-lives, probabilities of /^-delayed neu- 
tron emission Pn, neutron-capture cross sections, and ground-state spins and deformations are 
of importance in r-process calculations. ^^'^^) For the first time, most of these quantities can 
be obtained from a single, unified model, namely our present macroscopic-microscopic model. 
Neutron-capture rates are so far calculated with the statistical Hauser-Feshbach model, as long 
as the level density in the compound nucleus is sufficiently high to justify such an approach. 
As is the case for proton captures discussed in Sec. |5.1| , this method requires a knowledge of 
additional quantities, such as optical potentials and giant-dipole-resonance and level-density pa- 
rameters. For nuclei in the r-process path, in particular near closed neutron shells, the level 
density is small and the Hauser-Feshbach approach might no longer be applicable. In these cases, 
Breit-Wigner resonance capture and direct capture have to be considered, again requiring addi- 
tional nuclear-physics properties as input for the respective model calculations. ^^'^■^) For nuclei 
with Z > 80, fission barriers and rates of /3-delayed as well as neutron-induced fission are also 
important. Finally, within the recently favored "hot-entropy-bubble" r-process scenario, ^^'^^) 
charged-particle reactions during the so-called a-rich freezeout are also required in the ^ ~ 80 
region. Taken together, all these nuclear data for thousands of mostly unknown isotopes require 
a huge reaction network for "complete" r-process calculations. 

Primarily in order to facilitate these complicated and time-consuming calculations, since 
1957 many attempts to predict the N^^q distribution were based on the simplified assumption 
of the (n,7) ^ (Tjh) equilibrium concept. ^^'^^'^^) When assuming in addition a steady-fiow 
equilibrium of (5 decays, the prediction of r-abundances requires only the input of nuclear masses 
and corresponding neutron separation energies S'ln, /3-decay half-lives Tp and /3-delayed neutron- 
emission probabilities Pa, as well as the stellar parameters Tg, rin, and the process duration. 
Whereas for a given rin the 5in determine the r-process path, the Tf^ of the isotopes along this 
fiow path determine, in principle, the progenitor abundances and, when P^. branching during 
freezeout is considered, also the final r-abundances. Only in recent years could the validity of 
this "waiting-point" approximation in combination with a steady /3-decay flow be conflrmed 
locally for the >1 ~ 80 and 130 N^-^q peaks on the basis of the first experimental information in 
the r-process path. These recent results showed clearly, for example, how the long Tp of the 
classical = 82 waiting-point nucleus "'^^''Cd directly correlated with the large N^^q value of its 
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isobar ^^''Te in the A ~ 130 abundance peak seen in Fig. [T^ . 

Since there is not yet complete consensus on the stellar site of the r-process and the specific 
astrophysical conditions under which it takes place, one deductive approach to theoretical r- 
process studies has been to take the A^r,© observables as a constraint that allow the derivation 
of the necessary conditions required to reproduce these features. Using the unified nuclear- 
physics basis presented in this paper, supplemented by all experimental data available up to 
1991 as well as local improvements of the QRPA calculations, Kratz et al.^^) have obtained 
significant progress in reproducing the N^-^q pattern relative to the situation five or ten years 
ago. As an example, we show in Fig. 42 results of r-process calculations for two different mass 
models. The top calculation is based on our 1991 version of the FRDM.'^^'S'^) Our current 
FRDM (1992) shown in Fig. 19 has smaller mass and errors, but our preliminary r-process 
studies with it do not show significantly improved agreement. In the lower part of Fig. 
nuclear masses and input for the QRPA calculations were taken from the ETFSI-1 (1992) model 
in its preliminary, privately circulated version. 

With the current agreement between A^r,© data and the calculations, the theoretical treat- 
ment is sufficiently accurate that some conclusions can be drawn about the stellar conditions 
responsible for the production of r-process nuclei. For example, one has found that it is not 
possible to reproduce the A^r,© curve assuming a global steady-flow process. Instead, a mini- 
mum of three r-components with different neutron densities is required. Each of the components 
proceeds up to one of the abundance peaks and reaches a local steady-flow equilibrium which 
breaks down at the top of each peak, situated at one of the A'" = 50, 82, and 126 magic shells. 
With these results, the explosive He-burning scenarios favored in the 1980s can definitely be 
ruled out as possible sites for the r-process. 

The large deviations from the A^r,© pattern evident in Fig. especially those just before 
the abundance peaks at ^ ~ 130 and 195, were interpreted ^^'^^'®^) as arising from overly strong 
magic-neutron shell corrections as one moves away from a doubly magic configuration, from the 
neglect of the proton-neutron residual interaction, and from correlated problems with describing 
shape transitions in the neutron mid-shell regions around A^ = 66 and 104. 

Rather than representing a failure of the FRDM (1992) far from stability, these difficulties are 
due to normal inaccuracies that occur anywhere in the chart of the nuclides, both near stability 
and far from stability. They must be expected and considered normal in a model based on such 
a simple effective interaction as the one-body single-particle potential with a simple pairing 
residual interaction. Despite some deficiencies, the current models, when used appropriately, 
have been sufficiently accurate to considerably advance our understanding of several astrophysical 
processes and to identify specific nuclear-structure features of nuclei far from stability near magic 
shells and close to the neutron drip lines. These nuclei are normally inaccessible to experiment, 
but unique signatures of their nuclear structure are well preserved as differences in the observed 
A^r,© pattern relative to the pattern calculated with the FRDM (1992). Our identification of 
specific nuclear-structure features in nuclei far from stability has generated substantial interest in 
the nuclear-physics community and stimulated several calculations of nuclear masses in limited 
regions far from stability in terms of self-consistent mean-field theories. ™~^^) Some of these 
local calculations now produce nuclear masses that remove some of the discrepancies that are 



present in Fig. 42 



However, to be considered an improvement over current global nuclear-structure calculations, 
a new theory must achieve more than just a better local description of features that were 
postulated on the basis of existing theories. To justify the designation new and improved, a new 
theory should globally reproduce better nuclear masses and other ground-state nuclear-structure 
effects than does the FRDM (1992). It should achieve these results with a simple, global choice 
of constants. And, it should, just as has the current FRDM (1992), correctly predict new 
nuclear-structure features far from stability and serve to further enhance our understanding of 



p. MoUer, J. R. Nix, and K.-L. Kratz/Nuclear Properties 



27 



several astrophysical processes. The results achieved so far in new mean- field calculations are 
important preliminary steps in this direction. 

We have seen that certain differences in the abundance fits and remaining deficiencies can 
be attributed to the nuclear mass models applied. We now want to check whether effects from 
different half-life sets can also be separated out in our Nr^Q calculations. 

As we have discussed earlier, once the r-process path is defined by a contour line of 
constant 5in values, the abundances of isotopes in this path are directly related to their /3-decay 
half-lives according to the waiting-point concept. These /3-decay half-lives are determined by 
the energy window Q^, on the one hand, and by the low- lying nuclear structures in the Gamow- 
Teller strength function, on the other hand. As has been discussed in detail elsewhere, 
is predominantly influenced by the latter properties. Therefore, statistical models such as the 
gross theory of j3 decay, which neglect nuclear structures, will not be able to describe the 
/3-decay quantities far from stability in an adequate way. It has been well known for more than a 
decade that, due to the missing nuclear-structure effects, the /3-decay half-lives from this model 
are systematically too long by factors of five to 10 far from stability. Nevertheless, because the 
gross-theory predictions are still occasionally used in astrophysical calculations, ^4,74^ have for 
the purpose of illustrating its deficiencies performed some and Pn calculations with this model 
using Qp values from both the FRDM (1992) and ETFSI-1 (1992) model. We can show that the 
differences in the Qf^ predictions of the two models for isotopes in the regions of the r-process 
path affect the Tp values obtained in the gross theory only to a small extent. Therefore the 
effects in the calculated r-abundances originate mainly from the missing nuclear structures in 
the /3-strength functions of the gross theory. 

In Fig. ^ we present a comparison of r-process calculations based on the gross theory and 
on the QRPA. The upper part of Fig. ^ shows that for our best-fit conditions of the third 
component, ^^) a switch from the T^(QRPA) consistent with the FRDM masses to T^(g.t.) 
creates r-overabundances in the range 150 < A < 180. It is possible to adjust for this artificial 
effect by a change of astrophysical conditions, for example, by increasing in order to obtain a 
more neutron-rich r-process path and shorter T^(g.t.) in that mass region. ''^) However, as can 
be seen in the lower left-hand part of Fig. higher neutron densities, which improve the A^r,© 
fits in the 150 < A< 180 range, but are difficult to obtain in realistic r-process scenarios, ^^'^^) 
shift the structure of the A ~ 195 peak too far to lower masses and make such an approach 
invalid. We can only conclude that far from (3 stability the values of r^(g.t.) are, in fact, too 
long. 



6 SINGLE-PARTICLE LEVELS 

As a final result we present in Figs. 44-^ calculated proton and neutron single-particle level 
diagrams for representative spherical and deformed nuclei throughout the periodic system. The 
diagrams are useful for obtaining the spin and parity of low-lying states of odd-even nuclei 
and for identifying gaps in the level spectra that may be associated with particularly stable 
proton- neutron combinations. 

For heavier nuclei higher-multipole deformations become increasingly important, as can be 
seen in the mass table and color overview figures in our mass calculations.^). For prolate 
shapes we have therefore chosen €4 and eg to be functions of 62- The hexadecapole deformation 
€4 is either a linear function of 62 in the entire range < 62 < 0.4 or linear in the range 
< £2 < 0.2 and constant in the range 62 > 0.2, as indicated on the upper horizontal axis. 
The hexacontratetrapole deformation eg is zero when it is not explicitly mentioned in the figure 
captions. Otherwise, it varies linearly in the entire range < e2 < 0.4 in the manner indicated 
in the figure captions. 

One notes that proton number Z = 14 and neutron number A*" = 14 correspond to well- 
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developed spherical gaps for several of the lighter systems. The effect of the higher-multipole 
deformation parameters 64 and eg is clearly visible in the level-diagram sequences shown in 
Figs. ^^53 and 38-93, for example. 

In the heavy region it is particularly interesting that the deformed shell gaps at proton 
numbers Z = 104-110 and neutron number N = 162 emerge only for relatively large positive 
values of the hexadecapole deformation parameter €4. These gaps give rise to unusual stability, 
which has made possible the discovery of several new elements in this region. 26-30^ 

Figure ^ shows the predicted large spherical neutron gap = 184. However, the unusually 
large density of single-particle levels just above has the consequence that the largest negative 
microscopic correction occurs approximately at neutron number N = 178 in our model. 
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Figure captions 

Fig. |l] Global microscopic pairing gap for protons. The single-particle levels entering the 

pairing calculation correspond to ground-state shapes that have been determined 
by minimizing the total potential energy with respect to £2, £4, £3, and e% shape 
degrees of freedom. The jagged black lines indicate regions where experimental 
proton pairing gaps may be extracted from fourth-order odd-even mass differ- 
ences. Magic proton and neutron numbers are indicated by pairs of thin, parallel 
lines. In the Lipkin-Nogami model it is the sum Ap + A2p plotted here that should 
be compared to odd-even mass differences. In contrast to the behavior of BCS 
solutions, this sum shows no sign of collapse at magic proton numbers. 

Fig. |2| Global microscopic pairing gap for neutrons. The single-particle levels entering 

the pairing calculation correspond to ground-state shapes that have been deter- 
mined by minimizing the total potential energy with respect to £2, £4, £3, and 
£6 shape degrees of freedom. The jagged black lines indicate regions where ex- 
perimental neutron pairing gaps may be extracted from fourth-order odd-even 
mass differences. Magic proton and neutron numbers are indicated by pairs of 
thin, parallel lines. In the Lipkin-Nogami model it is the sum + A2n plotted 
here that should be compared to odd-even mass differences. In contrast to the 
behavior of BCS solutions, this sum shows no sign of collapse at magic neutron 
numbers. 

Fig. m Discrepancy between experimental proton pairing gaps determined from fourth- 

order odd-even mass differences and microscopic pairing gaps AlNp = Ap + A2p 
obtained in the Lipkin-Nogami model. Magic proton and neutron numbers are 
indicated by pairs of thin, parallel lines. Fairly large discrepancies occur in several 
places. However, these discrepancies do not necessarily mean that the calculated 
gaps are incorrect. Instead, it may be that the pairing gap is not determined 
properly from odd-even mass differences, as discussed in the text. In particular, 
when large errors in this figure occur in the same region where large errors occur 
in Fig. ^, sudden shape transitions are probably responsible. 

Fig. ^ Discrepancy between experimental neutron pairing gaps determined from fourth- 

order odd-even mass differences and microscopic pairing gaps AlNh = + A2n 
obtained in the Lipkin-Nogami model. Magic proton and neutron numbers are 
indicated by pairs of thin, parallel lines. Fairly large discrepancies occur in several 
places. However, these discrepancies do not necessarily mean that the calculated 
gaps are incorrect. Instead, it may be that the pairing gap is not determined 
properly from odd-even mass differences, as discussed in the text. In particular, 
when large errors in this figure occur in the same region where large errors occur 
in Fig. ^, sudden shape transitions are probably responsible. 

Fig. |5| Neutron separation energy for odd-neutron nuclei. Each odd-neutron nucleus 

is represented by a color field one unit high and two units wide. Black squares 
denote /^-stable nuclei. Where available, experimental masses were used to deter- 
mine the location of /3-stable nuclei; otherwise, calculated masses were used. The 
region of known nuclei is enclosed by a jagged black line, and magic proton and 
neutron numbers are indicated by pairs of thin, parallel lines. The region of the 
nuclear chart where the r-process occurs is approximately 2.0 MeV < 5in < 3.0 
MeV for odd A^, which in this figure appears in light blue. Recent experiments 
have reached the r-process region at both = 50 and iV = 82. 
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Neutron separation energy for even-neutron nuclei. Each even-neutron nucleus 
is represented by a color field one unit high and two units wide. Black squares 
denote /3-stable nuclei. Where available, experimental masses were used to de- 
termine the location of /3-stable nuclei; otherwise, calculated masses were used. 
The region of known nuclei is enclosed by a jagged black line, and magic proton 
and neutron numbers are indicated by pairs of thin, parallel lines. 

Two-neutron separation energy for odd- and even-neutron nuclei. Each nucleus is 
represented by a color field one unit high and one unit wide. Black squares denote 
/3-stable nuclei. Where available, experimental masses were used to determine the 
location of /3-stable nuclei; otherwise, calculated masses were used. The region 
of known nuclei is enclosed by a jagged black line, and magic proton and neutron 
numbers are indicated by pairs of thin, parallel lines. 

Discrepancy between experimental and calculated one-neutron separation ener- 
gies. Magic proton and neutron numbers are indicated by pairs of thin, parallel 
lines. Large discrepancies occur where there are large changes in the mass error 
between neighboring nuclei, namely at magic numbers and in the light region of 
nuclei. In the deformed rare-earth and actinide regions the discrepancy is very 
small. 

Proton separation energy for odd-proton nuclei. Each odd-proton nucleus is rep- 
resented by a color field one unit wide and two units high. Black squares denote 
/3-stable nuclei. Where available, experimental masses were used to determine 
the location of /3-stable nuclei; otherwise, calculated masses were used. The re- 
gion of known nuclei, which is enclosed by a jagged black line, extends in several 
places to where proton emission is energetically allowed from odd-proton nuclei. 
Magic proton and neutron numbers are indicated by pairs of thin, parallel lines. 

Proton separation energy for even-proton nuclei. Each even-proton nucleus is 
represented by a color field one unit wide and two units high. Black squares 
denote /^-stable nuclei. Where available, experimental masses were used to de- 
termine the location of /5-stable nuclei; otherwise, calculated masses were used. 
The region of known nuclei, which is enclosed by a jagged black line, nowhere ex- 
tends to where proton emission is energetically allowed from even-proton nuclei. 
Magic proton and neutron numbers are indicated by pairs of thin, parallel lines. 

Two-proton separation energy for odd- and even-proton nuclei. Each nucleus is 
represented by a color field one unit wide and one unit high. Black squares denote 
/5-stable nuclei. Where available, experimental masses were used to determine the 
location of /3-stable nuclei; otherwise, calculated masses were used. The region of 
known nuclei, which is enclosed by a jagged black line, barely extends to where 
two-proton emission is energetically allowed. Magic proton and neutron numbers 
are indicated by pairs of thin, parallel lines. 

Discrepancy between experimental and calculated one-proton separation ener- 
gies. Magic proton and neutron numbers are indicated by pairs of thin, parallel 
lines. Large discrepancies occur where there are large changes in the mass error 
between neighboring nuclei, namely at magic numbers and in the light region of 
nuclei. In the deformed rare-earth and actinide regions the discrepancy is very 
small. 
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Energy released in a decay. Black squares denote /3-stable nuclei. Where avail- 
able, experimental masses were used to determine the location of /3-stable nuclei; 
otherwise, calculated masses were used. The region of known nuclei is enclosed 
by a jagged black line, and magic proton and neutron numbers are indicated by 
pairs of thin, parallel lines. The calculations show that Qa is in the range 6-9 
MeV in the heaviest known region and 9-12 MeV in the deformed superheavy 
island surrounding ^^^110. 

Global a-decay half-life calculated from a semi-empirical relationship between 
Tq, and Qct- Black squares denote /5-stable nuclei. Where available, experimental 
masses were used to determine the location of /3-stable nuclei; otherwise, cal- 
culated masses were used. The region of known nuclei is enclosed by a jagged 
black line, and magic proton and neutron numbers are indicated by pairs of thin, 
parallel lines. 

Global microscopic /3-decay half- life for allowed Gamow- Teller transitions. Black 
squares denote /3-stable nuclei. Where available, experimental masses were used 
to determine the location of /3-stable nuclei; otherwise, calculated masses were 
used. The region of known nuclei is enclosed by a jagged black line, and magic 
proton and neutron numbers are indicated by pairs of thin, parallel lines. Above 
the black squares the calculated combined half-life with respect to (5'^ decay and 
electron capture is plotted; below the black squares the half-life with respect to 
/3~ decay is plotted. 

Features of the r-process. Black squares denote /3-stable nuclei. Where avail- 
able, experimental masses were used to determine the location of /^-stable nuclei; 
otherwise, calculated masses were used. The colored region in the main graph 
shows calculated half-life with respect to f3~ decay. The jagged black line gives 
the right-hand boundary of the region of known nuclei. The thick magenta line 
represents Si^ = 2.4 MeV, which is the approximate location of the r-process 
path for a particular set of stellar conditions. The magenta squares in the region 
of /3-stable nuclei are created in decay from the r-process line. The solar r-process 
abundance shown in the insert is plotted versus the mass number A, whose axis 
is curved slightly to follow the line of /3-stability. A line perpendicular to the 
valley of /3-stability and originating at a particular mass value crosses the A axis 
of the insert plot at right angles at this value and also passes through the circle 
giving the abundance for this A value. 

Comparison of calculated and experimental nuclear ground-state spins and pari- 
ties for odd-even nuclei in the light and medium-mass regions. Spherical assign- 
ments are used in the calculations when |e2| < 0.15. Many of the discrepancies 
occur in transition regions between spherical and deformed nuclei or where several 
levels are grouped close together. 

Similar to Fig. |l^ but for the heavy region. The discrepancies in the heaviest part 
of the actinide region occur because here several neutron single-particle levels are 
grouped very close together. 

Comparison of experimental and calculated microscopic corrections for 1654 nu- 
clei in the 1992 version of the finite-range droplet model. ^'^^) The bottom part 
showing the difference between these two quantities is equivalent to the difference 
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between measured and calculated ground-state masses. There are almost no sys- 
tematic errors remaining for nuclei with > 65, for which region the theoretical 
error is only 0.448 MeV. 

Calculated Gamow- Teller /3-strength function for ^^Rb. Whereas for the calcu- 
lation of the quantities in the Table we use the shapes obtained in our mass 
calculation, ^) in this figure we use the spherical shape because we wish to show 
a typical spherical /3-strength function and because experimentally ^^Rb is spher- 
ical. The arrows with wide heads denote successive neutron-separation energies 
in the daughter, and the arrow with a thin head denotes the value oi Qp. The 
calculated strength containing only a few large peaks in the low-energy region is 
typical of a spherical nucleus. Since all the strength is calculated to be above the 
one-neutron separation energy, the theoretical probability P^ for /3-delayed neu- 
tron emission is 100%. Experimentally, there is a large peak in the GT strength 
function in the region 3.5-4.0 MeV, so experimentally the /3-delayed neutron 
emission probability is only 8.5%. 

Calculated Gamow- Teller /J-strength function for ^^Rb. Here we use the shape 
obtained in our mass calculation. The strength function is typical of that of a 
deformed nucleus. Because there is a higher likelihood of significant strength 
in the low-energy region for deformed nuclei than for spherical nuclei there is 
a characteristic, large decrease in the /?-decay half-life at the shape transition. 
Since there is considerable strength below the neutron-emission threshold the 
/3-delayed neutron-emission probability is low. 

Ratios between calculated and experimental half-lives for /5~ decay as functions 
of neutron number N . There are no systematic effects versus N . 

Ratios between calculated and experimental half-lives for /3~ decay as functions 
of the experimental half-life for f5~ decay. As expected, we find a very strong 
correlation between the error and the experimental /3-decay half-life. The cor- 
relation is such that we can expect fairly reliable half-life calculations far from 
/3-stability, in the region of interest for astrophysical r-process calculations. An 
analysis of the results in this figure is presented in Table |2| and is also discussed 
in the text. 

Ratios between calculated and experimental half-lives for /5~ decay as functions 
of Qp. The discrepancy is expected to be larger for low values of Qp because the 
calculated half-life is most sensitive here to errors in the positions of the peaks 
in the strength functions. 

Ratios between calculated and experimental half-lives for I3~^ decay and electron 
capture as functions of neutron number N . There are no systematic effects versus 
N. 

Ratios between calculated and experimental half-lives for /?"*" decay and electron 
capture as functions of the experimental half-life for /?"'" decay and electron cap- 
ture. Surprisingly, there is no strong correlation between the error and TJg^exp; 
except for odd-odd nuclei. An analysis of the results in this figure is presented 
in Table ^ and is also discussed in the text. 

Ratios between calculated and experimental half-lives for /J"*" decay and electron 
capture as functions Qp. Only the error for odd-odd nuclei is correlated with 
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Comparison between experimental and calculated energy releases Qa for the 
N = 154 and 155 isotonic chains. The experimental data are from Fig. 6 of 
Miinzenberg et al.,^'^) where the data are compared to predictions of the FRLDM 
(1988). The new results presented here agree with the data much better. The 
improvement is due partly to the inclusion of the eg shape degree of freedom. 

Calculation to show the reliability of the FRDM (1992) in new regions of nuclei. 
The FRDM (1992) was adjusted to 1654 masses known in 1989. 2) The figure 
shows the deviations between experimental and calculated masses for 217 new 
nuclei whose masses were measured between 1989 and 1993.^'^) The error is 
4% smaller in the new region compared to that in the region where the model 
constants were adjusted. There are no systematic effects visible in the figure. 

Calculation to show the reliability of the ETFSI-1 (1992) in new regions of nuclei. 
The ETFSI-1 (1992) was adjusted to masses known in 1988.^^) The fi gure shows 
the deviations between experimental and calculated masses for 210 new nuclei 
whose masses were measured between 1989 and 1993. ^'^) The error is 10% larger 
in the new region compared to that for the 1989 data set^) we have available, 
which is only marginally different from the 1988 data set to which the model 
constants were adjusted. There are no systematic effects visible in the figure. 

Similar to Figs. ^ and ^ but for the mass formula of von Groote et al. ^°) With 
its postulated shell corrections and more adjustable constants than in our model 
with calculated shell corrections, the error in the new region is 104% larger than 
for a 1977 set of measured masses, ^^) which is only marginally different from 
the set of masses where the constants were adjusted. There is also a systematic 
increase in the error with increasing distance from /? stability. 

Calculation to show the reliability of one-neutron separation energies obtained 
from the FRDM (1992) in new regions of nuclei. There are no systematic effects 
visible in the figure. 

Calculation to show the reliability of one-neutron separation energies obtained 
from the ETFSI-1 (1992) model in new regions of nuclei. There are no systematic 
effects visible in the figure, apart from an odd-even staggering related to problems 
in the pairing part of the ETFSI-1 (1992) model. 

Calculation to show the reliability of one-neutron separation energies obtained 
from the mass formula of von Groote et al. in new regions of nuclei. On the 
neutron-rich side the calculated one-neutron separation energies are systemati- 
cally too high by 0.29 MeV. 

Calculation to show the reliability of energy releases for /3~ decay obtained from 
the FRDM (1992) in new regions of nuclei. There are no systematic effects visible 
in the figure. 

Calculation to show the reliability of energy releases for (3~ decay obtained from 
the ETFSI-1 (1992) model in new regions of nuclei. 

Calculation to show the reliability of energy releases for f3~ decay obtained from 
the mass formula of von Groote et al. in new regions of nuclei. On the neutron- 
rich side the calculated energy releases are systematically too low by 0.52 MeV. 
Corresponding /?-decay rates based on these values would be too slow. 
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Calculation to show the reliability of the TF (1994) model in new regions of 
nuclei. The TF (1994) model was adjusted to 1654 masses known in 1989.^) 
The figure shows the deviations between experimental and calculated masses for 
217 new nuclei whose masses were measured between 1989 and 1993.^'"^) The 
error is 15% larger in the new region compared to that in the region where the 
model constants were adjusted. There are no systematic effects visible in the 
figure. 

Calculation to show the reliability of the Wapstra-Audi systematic-trend mass 
model in new regions of nuclei. The systematic-trend masses were provided in 
the 1989 midstream mass evaluation. ^) The figure shows the deviations between 
experimental and systematic masses for 187 new nuclei that were given by sys- 
tematic trends in the 1989 evaluation and whose masses were measured between 
1989 and 1993.^'^) There are no systematic effects visible in the figure. 

Comparison between energy releases Qa obtained in the FRDM (1992), ETFSl-1 
(1992) model, and FDSM (1992) and recent experimental data for the heaviest 
known element. '^^) When several values of Qa were measured we choose for the 
figure the highest value. 

The hot hydrogen burning cycles, ^'') typically consisting of three proton captures, 
two /3+ decays, and a closing (p,a) reaction. Break-out towards heavier nuclei 
occurs only via the (p,7)/(p,a) branching at the cycle closings. Because of the 
^^F(p,a) reaction no OFNe cycle exists, which would otherwise connect the CNO 
and NeNaMg cycles. 

Calculated r-process abundances (solid lines) compared to measured values (solid 
circles). For both the upper and lower parts of the figure /3-decay half-lives and 
delayed-neutron emission probabilities are calculated in a QRPA model based on 
folded- Yukawa single-particle energies, but experimental information has been 
used when available. In the upper part of the figure the r-process path was 
determined from the FRDM (1991),^^'^^) and in the lower part of the figure it 
was determined from the preliminary, privately circulated version of the ETFSl-1 
(1992) model. ^^'^8) 

Static steady-flow calculations of the r-process abundance N^-^q for the 135 < 
A < 195 mass region, the so-called "third-component." The right-hand part of 
the figure shows the best fit obtained with and Pn values from our QRPA 
model. The upper left-hand part of the figure shows N^^q for the same Tg-rin 
conditions obtained with the on-the-average five-times longer values from the 
gross theory. ''^) As can be seen from the lower left-hand part of the figure, with 
these Tjs values reasonable flts for the actinide region require neutron densities 
of 10^^ cm~^, which are difficult to obtain in the hot-entropy-bubble r-process 
scenario. 

Calculated proton single-particle level diagram for nuclei in the vicinity of ^gOg. 
Calculated neutron single-particle level diagram for nuclei in the vicinity of ^gOg. 



Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 
ity of 2oCai4. 



Calculated neutron single-particle level diagram for proton-rich nuclei in the 
vicinity of 2oCa;^4. 
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Calculated proton single-particle level diagram for nuclei in the vicinity of 2oCa24. 
Calculated neutron single-particle level diagram for nuclei in the vicinity of 

44p 

Calculated proton single-particle level diagram for neutron-rich nuclei in the 
vicinity of 2t^^M- 

Calculated neutron single-particle level diagram for neutron-rich nuclei in the 
vicinity of 2oCa34. 

Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 
ityof iNi22. 

Calculated neutron single-particle level diagram for proton-rich nuclei in the 
vicinity of 28^122- 

Calculated proton single-particle level diagram for nuclei in the vicinity of 28^136- 
Calculated neutron single-particle level diagram for nuclei in the vicinity of 

28^^136- 



Calculated proton single-particle level diagram for neutron-rich nuclei in the 
vicinity of ssNigo- 

Calculated neutron single-particle level diagram for neutron-rich nuclei in the 
iNiso- 



vicinity of 28^15 



Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 
ity of 4oZr4o. The hexacontratetrapole deformation eg is zero for oblate shapes 
and is given by eg = — 0.1e2 for prolate shapes. 

Calculated neutron single-particle level diagram for proton-rich nuclei in the 
vicinity of 4QZr4Q. The hexacontratetrapole deformation e% is zero for oblate 
shapes and is given by eg = — 0.1e2 for prolate shapes. 

Calculated proton single-particle level diagram for nuclei in the vicinity of 4oZr5o. 
Calculated neutron single-particle level diagram for nuclei in the vicinity of 

90 7-, 
40 ■^^50- 

Calculated proton single-particle level diagram for neutron-rich nuclei in the 
vicinity of ^4QZrgg. The hexacontratetrapole deformation eg is zero for oblate 
shapes and is given by eg = — 0.1e2 for prolate shapes. 

Calculated neutron single-particle level diagram for neutron-rich nuclei in the 
vicinity of ^4oZrgg. The hexacontratetrapole deformation eg is zero for oblate 
shapes and is given by eg = — 0.1e2 for prolate shapes. 



Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 
ity of loosnso. 

Calculated neutron single-particle level diagram for proton-rich nuclei in the 

00c„ 
50^^50- 



vicinity of ^goSugQ. 



Calculated proton single-particle level diagram for nuclei in the vicinity of ^5oSngg. 
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Fig. ^ Calculated neutron single-particle level diagram for nuclei in the vicinity of 

116c„ 

5oongg. 



Fig. ^ Calculated proton single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^ioSnga- 

Fig. ^ Calculated neutron single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^5oSn82. 



Fig. ^ Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 

ity of ^iSmyg. 



Fig. ^ Calculated neutron single-particle level diagram for proton-rich nuclei in the 

vicinity of ^02811170 . 

Fig. ^ Calculated proton single-particle level diagram for neutron-rich nuclei in the 

vicinity of "^glSnige- 

Fig. ^ Calculated neutron single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^g^Smgg. 

Fig. ^ Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 

ity of if^Wg^. 



Fig. ^ Calculated neutron single-particle level diagram for proton-rich nuclei in the 

^82- 



vicinity of ^74W8 



Fig. ^ Calculated proton single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^74W]^;^g. 

Fig. ^ Calculated neutron single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^y^W^^^g. 



Fig. ^ Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 

ity of liPbgg. 

Fig. ^ Calculated neutron single-particle level diagram for proton-rich nuclei in the 

vicinity of ^I^Pbgg. 

Fig. ^ Calculated proton single-particle level diagram for nuclei in the vicinity of ^82^^126 • 



Fig. ^ Calculated neutron single-particle level diagram for nuclei in the vicinity of 

208pK 
82-'^"l26- 

Fig. ^ Calculated proton single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^giPb^sg, which is located on the r-process path and is calculated 
to be deformed in its ground state. 

Fig. ^ Calculated neutron single-particle level diagram for neutron-rich nuclei in the 

vicinity of ^gifbisg, which is located on the r-process path and is calculated to 
be deformed in its ground state. 



Fig. ^ Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 

ity of 2i6Thi2g. 

Fig. ^ Calculated neutron single-particle level diagram for proton-rich nuclei in the 

90 -'^^^126 • 



vicinity of ^ggTh^^ 
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Calculated proton single-particle level diagram for nuclei in the vicinity of ^gQTh^42. 
The hexacontratetrapole deformation eg is zero for oblate shapes and is given by 
eg = 0.1e2 for prolate shapes. 

Calculated neutron single-particle level diagram for nuclei in the vicinity of 
^gQTh^42. The hexacontratetrapole deformation ee is zero for oblate shapes and 
is given by eg = 0.1 e2 for prolate shapes. 

Calculated proton single-particle level diagram for proton-rich nuclei in the vicin- 
ity of fQQFm^4Q. The hexacontratetrapole deformation eg is zero for oblate shapes 
and is given by eg = 0.1 5e2 for prolate shapes. 

Calculated neutron single-particle level diagram for proton-rich nuclei in the 
vicinity of fooFiiii4o- hexacontratetrapole deformation eg is zero for oblate 
shapes and is given by eg = 0.15e2 for prolate shapes. 

Calculated proton single-particle level diagram for nuclei in the vicinity of f ooFm^52 
The hexacontratetrapole deformation eg is zero for oblate shapes and is given by 
eg = 0.2e2 for prolate shapes. 

Calculated neutron single-particle level diagram for nuclei in the vicinity of 
^QQFm^52- The hexacontratetrapole deformation eg is zero for oblate shapes and 
is given by eg = 0.2e2 for prolate shapes. 

Calculated proton single-particle level diagram for nuclei in the vicinity of f QoFm^g2 
Calculated neutron single-particle level diagram for nuclei in the vicinity of 

100-^™162- 

Calculated proton single-particle level diagram for superheavy nuclei in the vicin- 
ity of ^'^^110^52, which is calculated to be deformed in its ground state. 

Calculated neutron single-particle level diagram for superheavy nuclei in the 
vicinity of ^''■^110]^g2, which is calculated to be deformed in its ground state. 

Calculated proton single-particle level diagram for superheavy nuclei in the vicin- 
ity of ^^^114]^g4, which is calculated to be spherical in its ground state. 

Calculated neutron single-particle level diagram for superheavy nuclei in the 
vicinity of ^^^114]^g4, which is calculated to be spherical in its ground state. 



Calculated proton single-particle level diagram for superheavy nuclei in the vicin- 

184) 



ity of '^'^^124^g4, which is calculated to be spherical in its ground state. 



Calculated neutron single-particle level diagram for superheavy nuclei in the 
vicinity of ^^^12^^^^, which is calculated to be spherical in its ground state. 
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Table. Calculated Nuclear Ground-State Properties 



Z Proton number. The Table is ordered by increasing proton number. The corre- 

sponding chemical symbol of each named element is given in parentheses. For 
consistency we use the same naming scheme as in our earlier publication,^) al- 
though the names of some of the heavier elements may be officially changed in 
the future. 

N Neutron number. 

A Mass number. 

Qp Projection of the odd-proton angular momentum along the symmetry axis and 

parity of the wave function. 

OJJ Projection of the odd-neutron angular momentum along the symmetry axis and 

parity of the wave function. 

AlNp Pairing gap for protons in the Lipkin-Nogami model, given by Ap + A2p. 

AlNi, Pairing gap for neutrons in the Lipkin-Nogami model, given by An + A2n- 

-E'bind Total binding energy. 

Sill Onc-ncutron separation energy. 

S211 Two-neutron separation energy. 

Pa Probability for producing a final nucleus with mass number A following /? decay 
and delayed neutron emission. 

Pa-1 Probability for producing a final nucleus with mass number A — 1 following P 

decay and delayed neutron emission. 

Pa-2 Probability for producing a final nucleus with mass number A — 2 following P 

decay and delayed neutron emission. 

Qp Energy released in /3 decay. 

Tp Half-life with respect to Gamow-Teller /3 decay. 

Sip One-proton separation energy. 

S'2p Two-proton separation energy. 

Qa Energy released in a decay. 

Ta Half-life with respect to a decay. 
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